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1.0  DYSCO  OVERVIEW 

DYSCO  Is  a  fully  interactive  computer  program  for  modeling  arbitrarily  coupled 
dynamic  and  aerodynamic  systems.  It  Is  designed  to  provide  the  automatic 
dynamic  coupling  of  components  and  associated  forces  necessary  to  formulate 
the  system  equations  of  motion  of  a  model  and  to  execute  user-selected  solu- . 
tion  algorithms.  The  program  is  composed  of  an  executive  which  acts  on  user- 
oriented  modeling  commands  to  control  execution,  a  data  base  management  facil¬ 
ity  for  Input/output  and  data  storage  and  retrieval,  and  an  expandable  tech¬ 
nology  library  of  component,  force,  and  solution  algorithms. 

A  component  technology  module  ("component")  is  a  group  of  program  subroutines 
that  represents  a  system  of  2nd-order  ordinary  differential  equations  which 
may  have  constant,  time-dependent,  or  nonlinear  coefficients.  The  coeffi¬ 
cients  may  be  arbitrary  (but  computable)  functions  of  the  independent  variable 
(time)  and  the  system  state  vector.  Thus,  a  "component"  is  an  algorithm  which 
defines  the  degrees  of  freedom  and  computes  the  coefficients  of  a  system  of 
differential  equations  from  specified  physical  parameters,  time,  and  the  com¬ 
ponent  state  vector.  A  force  technology  module  ("force")  is  an  algorithm 
which  computes  the  generalized  force  acting  on  the  component  as  a  function  of 
the  present  or  a  past  state  vector  of  the  component  and  other  physical  data 
(e.g.,  aerodynamic).  A  "model"  is  a  coupled  system  of  differential  equations 
formulated  from  the  component  equations  and  associated  forces  specified  by  the 
user.  A  solution  technology  module  can  include  any  algorithm  which  operates 
on  a  system  of  differential  equations  (e.g.,  eigenanalysis,  time  history). 

The  user  defines  a  model  by  selecting  components  and  forces  and  identifying 
the  data  to  be  used.  From  the  model  definition,  the  coupled  system  equations 
are  automatically  formed  independent  of  any  solution  algorithm.  Solutions  may 
then  be  selected  for  the  model.  Modeling  will  typically  be  accomplished  over 
several  Interactive  sessions  using  previous  results  as  building  blocks  for 
forming  variations  In  the  models.  Component  and  force  input  data  are  Individ¬ 
ually  stored  on  permanent  user  data  files.  An  existing  input  data  set  for  a 


1 


component  or  force  can  be  edited  to  form  a  different  data  set  or  to  replace 
the  original  data  set.  Model  definitions  are  also  Individually  stored  and  may 
be  edited  to  replace,  delete,  Insert,  or  add  components  and  forces  for  new 
model  formulation. 

DYSCO  supports  two  types  of  data  files:  random  access  user  data  files  and 
sequential  external  data  files.  A  user  data  file  contains  the  data  resulting 
from  Interactive  user  modeling  sessions  and  includes  component,  force,  and 
model  definitions  and  data,  and  auxiliary  engineering  data  such  as  airfoil 
tables.  The  data  on  user  data  files  is  in  a  format  that  allows  it  to  be 
recognized,  listed,  and  appropriately  manipulated  by  the  executive.  Examples 
of  external  data  are  induced  velocity  maps  for  rotor  systems,  airfoil  tables, 
and  plot  files.  While  the  executive  has  knowledge  of  and  tracks  the  existence 
of  external  data  during  a  session,  it  does  not  automatically  attempt  to  handle 
the  data.  Special  programming  must  be  developed  within  technology  modules  to 
utilize  the  data  as  required.  Some  data,  such  as  airfoil  tables,  although 
created  externally  and  initially  Introduced  as  an  external  data  file,  are  used 
frequently  and  are  thus  read,  reformatted,  and  saved  on  a  user  data  file  via  a 
special  modeling  command. 

1.1  MATHEMATICS  OF  THE  COMPONENT  AND  SYSTEM  EQUATIONS 

1.1.1  Assumptions.  The  general  concept  upon  which  DYSCO  is  based  depends  on 
the  following  assumptions: 

1.  The  relevant  physics  of  a  system  may  be  modeled  as  a  set  of 
second-order  differential  equations  in  the  time  domain.  The 
equations  are  of  the  form 

M  • 

MX  +  CX  +  KX  -  F 

where  X  is  the  vector  of  the  degrees  of  freedom  of  the  system. 

M,  C,  and  K  are  coefficient  matrices  and  F  is  a  force  vector 

where  M,  C,  K,  and  F  are  arbitrary  functions  of  X,  X,  and  time. 
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2.  It  Is  possible  to  formulate  the  equations  of  a  system  based  on 
the  equations  of  the  components  of  the  system.  The  equations 
of  the  components  are  of  the  same  form  as  the  equations  of  the 
system  as  in  1,  above.  "Formulation"  includes  the  establish¬ 
ment  of  the  logic  which  will  allow  the  computation  of  the  vary¬ 
ing  coefficients  (elements  of  M,  C,  and  K)  and  forces  (F)  of 
the  system  equations  based  on  the  computed  coefficients  and 
forces  of  each  component  at  any  point  in  time. 

3.  It  is  possible  to  compute  the  state  vector  (X  and  X)  of  each 
component  based  on  the  state  vector  of  the  system.  The  state 
vectors  of  the  components  may  be  used  to  compute  the  coeffi¬ 
cients  and  forces  of  the  equations  of  the  components  at  each 
point  in  time.  These  coefficients  and  forces  may  then  be  used 
to  compute  the  coefficients  and  forces  of  the  system  equations 
as  in  2,  above. 

1.1.2  Component  Equations.  Any  portion  of  the  system  which  can  be  repre¬ 
sented  by  a  set  of  second-order  differential  equations,  subject  to  the  condi¬ 
tions  stated  below,  may  be  considered  to  be  a  "component"  represented  in  its 
own  local  coordinate  system.  The  equation  of  component  I  may  be  written 

H,  +  C,  X,  +  K,  X,  -  F,  ♦  F,r 

where  Xj  is  a  vector  containing  the  degrees  of  freedom  of  the  component.  The 
coefficient  matrices,  Mj,  Cj,  and  Kj,  and  the  force  vector,  F j ,  may  be  arbi¬ 
trary  functions  of  Xj,  Xj,  and  time.  FjR  is  a  vector  of  the  reaction  forces  at 
the # interfaces  to  other  components.  The  elements  of  Fj  may  also  be  functions 
of  X  and  X  of  other  components.  "Arbitrary  function"  as  used  above  means  that 
it  is  possible  to  write  a  program  to  compute  the  numerical  value  of  the  func¬ 
tion  at  any  point  in  time  given  numerical  values  of  the  appropriate  variables. 
This  functionality  includes  linear,  nonlinear  analytical,  periodic,  and  tabu¬ 
lar  functions. 
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The  individual  degrees  of  freedom  must  represent  independent  (generalized) 
coordinates.  TfiS1' degrees  of  freedom  may  be  actual  physical  displacement*  in 
coordinate  systems  which  are  fixed  or  moving  at  a  uniform  velocity  or  rotating 
at  a  constant  velocity.  They  may  also  represent  modal  displacements  or  any 
other  generalized  coordinate.  Xj  may  contain  any  combination  of  the  allowable 
types  of  degrees  of  freedom. 

The  coupling  treated  in  DYSCO  involves  equating  physical  displacements  of  two 
or  more  components.  These  physical  displacements  may  or  may  not  be  general¬ 
ized  coordinates  of  the  components;  however,  they  must  be  expressible  as  a 
linear  combination  of  degrees  of  freedom  of  each  of  the  components  to  be  cou¬ 
pled.  Thus,  the  vector,  Xj,  must  contain  appropriate  degrees  of  freedom  to 
allow  desired  coupling  to  other  components. 

1.1.3  System  Equations.  The  system  coefficient  matrices  and  generalized 
force  vector  are  formed  from  the  component  matrices  and  vectors.  The  degree 
of  freedom  vector  of  the  system,  X,  must  contain  only  independent  coordinates. 
This  vector  is  formed  by  assembling  the  degrees  of  freedom  of  all  of  the  com¬ 
ponents  and  eliminating  duplications  when  specific  degrees  of  freedom  are 
coupled  or  eliminating  a  degree  of  freedom  whenever  linear  combinations  of 
degrees  of  freedom  are  coupled.  Also,  the  reaction  forces  between  components 
cancel  when  the  system  equations  are  assembled.  Thus,  they  may  be  neglected. 

1.1.4  Coordinate  Transformations.  Each  degree  of  freedom  of  a  component  is 
equal  to  a  degree  of  freedom  of  the  system  or  is  a  linear  combination  of  de¬ 
grees  of  freedom.  This  relationship  may  be  expressed  as  a  matrix  equation 


where  Tj  is  an  njxn  matrix,  nj  and  n  are  the  number  of  degrees  of  freedom  of 
the  component  and  the  system.  Since  the  coupling  treated  in  DYSCO  physically 
joins  components,  Tj  may  not  be  a  function  of  time  and,  therefore, 


When  the  system  contains  degrees  of  freedom  in  fixed  and  moving  coordinate 
systems,  any  necessary  transformations  are  to  be  modeled  within  the  component 
representations.  This  will  result  in  time  dependent  coefficients,  which  is 
consistent  with  the  component  equations  as  previously  defined. 


Using  the  transformation  relationships,  the  equation  for  a  component  can  be 
written  in  terms  of  the  System  degrees  of  freedom 


Hj  Tj  X  +  Cj  T,  X  +  Kj  Tj  X  ■  F, 

This  matrix  equation,  however,  is  not  a  proper  dynamic  representation  since  . 
the  coefficient  matrices  are  rectangular  and  Fj  is  not  the  generalized  force 
vector  corresponding  to  X.  The  geheralized  force  vector  is  defined  as  the 

vector  that  yields  work  done  by  the  external  forces  when  the  product  xl  FT  is 

T  1  A 

formed  where  Xj  is  Xj  transpose. 

The  work  done  cannot  depend  on  the  coordinate  system;  that  is,  the  work  due  to 
the  component  degrees  of  freedom  must  be  equal  to  the  work  due  to  the  system 
degrees  of  freedom.  Then,  since 


XT  p 
*1  hI 


XT  tt  f 
x  Ij  i-j 


the  right  side  of  the  transformed  matrix  equation  must  be  t|  Fj.  This  is 
achieved  by  premultiplying  by  t|,  yielding 


TI  MI  TI  X  +  TI  CI  TI  X  ♦  TI  KI  TI  x  ■  TI  FI 


When  this  transformation  is  performed  on  each  component  set  of  equations,  they 
are  consistent  and  may  be  added  to  form  the  equations  of  the  coupled  system 
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MX  +  cx  +  tgc  -  F 

where 

H  ■  pi  Hi  Ti 

c  -  ?  TI  ci  T> 

i  Ti 

f-ztJfi 

These  equations  are  not  limited  to  linear  equations.  It  should  be  noted,  how¬ 
ever,  that  when  the  coefficients  and  forces  are  functions  of  time  or  are  non¬ 
linear,  they  must  be  evaluated  numerically  at  discrete  values  of  time,  or  when 
the  component  state  vectors  have  been  evaluated  if  the  coefficients  or  forces 
are  nonlinear. 

(A  global  reference  system,  which  permits  the  application  of  consistent  gravi¬ 
tational  and  centrifugal  force  vectors  to  component  degrees  of  freedom,  has 
been  developed  as  an  adjunct  to  several  component  and  solution  technology  mod¬ 
ules  and  is  discussed  in  paragraph  3.4,  Volume  II.) 

1.2  SOLUTION  TECHNIQUES 

The  solutions  currently  available  in  DYSCO  are  of  three  types: 

1.  Eigenanalysis 

2.  Frequency  domain 

3.  Time  history. 

The  eigenanalysis  and  frequency  domain  solutions  are  purely  mathematical  oper¬ 
ations  on  the  system  constant  M,  C,  and  K  matrices.  Time  history  solutions 
involve  solutions  of  the  differential  equations  in  the  time  domain. 
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A  time  history  solution  is  an  integration  algorithm  which  takes  the  accelera¬ 
tions  at  a  given  time,  t,  and  computes  the  velocities  and  displacements  at  the 
next  time  increment,  t  +  At.  The  acceleration  is  obtained  by  solving 

X  -  M"1  (F  -  CX  -  KX) 

at  time,  t.  The  values  of  X  and  X  were  obtained  on  the  previous  execution  of 
the  integration  and  M,  C,.K,  and  F  are  obtained  as  previously  computed  con¬ 
stants  and  time-dependent  values,  if  any. 

The  time-dependent  component  coefficients  are  computed,  transformed,  and  added 
to  the  system  constant  matrices  at  each  time  step.  The  calculation  of  nonlin¬ 
ear  coefficients  and  forces  requires  information  derived  from  local  velocities 
and  displacements  of  the  component  (and  possibly  other  components).  This  in¬ 
formation  is  obtained  from  the  system  state  vector  by  use  of  the  transforma¬ 
tion  data  as  previously  described. 

1.3  TECHNOLOGY  MODULE? 

The  technology  module  library  contains  the  component,  force,  and  solution  al¬ 
gorithms  available  for  formulating  and  analyzing  a  model.  Each  technology 
module  is  a  group  of  program  subroutines  called  technical  modules  that  imple¬ 
ment  a  given  algorithm.  The  technology  module  library  is  designed  for  expan¬ 
sion,  allowing  new  methods  to  be  added  in  a  consistent  and  well-defined  manner. 
A  technology  module  has  multiple  Interfaces  with  the  executive  corresponding 
to  specific  steps  in  the  modeling  process  and  is-  thus  composed  of  a  specific 
set  of  technical  modules,  depending  on  whether  it  represents  a  component, 
force,  or  solution.  Technology  modules  are  given  a  four-character  name.  The 
first  character  is  C,  F,  or  S,  depending  on  whether  the  technology  module  is  a 
component,  force,  or  solution.  The  second  and  third  characters  are  used  as  a 
description  and  the  fourth  character  indicates  a  general  level  of  complexity. 
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Technical  Modules  -  The  type  and  number  of  technical  modules  required  for  1m- 
plementatlon  of  “a  technology  module  depend  on  the  technology  module  type. 
Each  technical  module  has  specific  functions  to  perform  although  the  data  re¬ 
quired,  the  algorithms,  and  the  output  data  vary.  The  name  of  a  technical 
module  Is  five  characters  -  the  technology  module  name  plus  a  character  that 
indicates  the  function  to  be  performed.  A  brief  description  of  the  required 
technical  modules  Is  shown  below  (full  detail  may  be  found  in  Volume  II  of 
this  report): 

Component  Technical  Modules 


Input  C---I 

Definition  C---D 

Coefficient  C---C 

Active  C — A 


Block 

C— B 

Loads 

C — L 

Defines  the  user  Input  data  requirements 

Defines  the  degrees  of  freedom  and  Implicit 
relationships  for  a  particular  usage 

Computes  the  constant  coefficients  and 
forces  of  the  component  equations 

Computes  time  varying  coefficients  and 
forces  of  the  component  equations  based  on 
component  state  vector.  Accesses  optional 
"force  module"  for  computation  of  applied 
forces. 

Reads/writes  private  common  block  /C---/ 

Optional.  Computes  time  history  component 
loads  based  on  component  state  vector. 


Force  Technical  Modules 


Input  F---I 
Coefficient  F---C 
Active  F---A 

Block  F---B 


Defines  the  user  Input  data  requirements 

Performs  preliminary  calculations 

Computes  time  varying  applied  forces  of  the 
associated  component  equation 

Reads/writes  private  common  block  /F — / 


\  . 

Solution  Technical  Modules 

Input  S---I  Defines  the  user  Input  data  requirements 

\ 

Active  S---A  Performs  the  solution  and  outputs  the 

results 

C---B  and  F---B  are  utilities  for  the  preparation  of  a  ’'case,"  which  is  auto¬ 
matically  formed  as  a  user  option.  A  case  includes  a  model,  user  input  for  a 
selected  solution,  and  any  required  auxiliary  data  (e.g.,  airfoil  tables). 
Thus,  the  CASE  command  allows  the  user  to  execute  a  solution  which  has  been 
previously  defined. 

During  an  interactive  session,  degree  of  freedom  names  are  either  automati¬ 
cally  formed  in  a  standard  format  or  they  may  be  supplied  by  the  user,  depend¬ 
ing  on  the  formulation  of  the  technology  module  for  a  given  component.  Com¬ 
ponents  are  automatically  coupled  Into  a  model  either  by  recognition  of  iden¬ 
tical  degree  of  freedom  names  or  by  generated  or  supplied  linear  constraint 
relations.  The  degrees  of  freedom  may  represent  physical  displacements  or 
generalized  modal  displacements.  ach  component  is  represented  independently, 
the  only  requirement  being  that  necessary  interface  coordinates.be  included  as 
degrees  of  freedom. 

1.4  COMMUNICATION 

Mechanisms  are  provided  for  communication  between  the  executive  and  technology 
modules,  between  technology  modules,  and  within  technology  modules. 

Between  the  Executive  and  Technology  Modules  -  The  executive  controls  data 
related  to  a  model  being  formed  or  executed  with  a  chosen  solution.  Technical 
modules  may  require  this  data  as  input  or  may  be  required  to  compute  it  and 
send  it  to  the  executive  for  storage.  These  requirements  are  specific  to  the 
particular  technical  module.  A  design  feature  which  prohibits  the  passage  of 
physical  data  into  the  executive  allows  the  independent  implementation  of 


9 


f 


technology  nodules.  All  physical  data  associated  with  a  technology  module  are 
treated  as  loca.1  data,  and  only  abstract  data  (e.g.,  equation  coefficients) 
are  passed  to  the  executive.  The  nechan Isns  utilized  are: 


1.  Argument  Lists 

2.  XBG  Utilities 


3.  /AFTAB/ 


-  Provide  Input  from  the  Executive  and  output 
to  the  Executive. 

-  Provide  for  retrieval  of  Base  and  Global 
Variables  (user  Input  for  Component  and 
Force  modules). 

-  A  common  block  used  by  the  Executive  to 
store  all  airfoil  tables  required  for  a 
model.  /AFTAB/  Is  Incorporated  Into  specif¬ 
ic  Technical  Modules  for  accessing  the 
tables. 


Between  Technology  Modules  •  Communication  between  technology  modules  Is  via 
special  "shared"  common  blocks.  Communication  between  the  technical  modules 
of  a  particular  component  or  between  a  component  and  a  force  or  solution  mod¬ 
ule  is  through  a  shared  common  block.  The  name  of  a  shared  common  block  is 
descriptive  of  its  usage. 


Within  a  Technology  Module  -  Communication  strictly  between  the  technical  mod¬ 
ules  for  a  specific  technology  module  is  through  a  "private"  common  block. 
The  name  of  a  private  common  block  is  the  same  as  the  component,  force,  or 
solution  module  with  which  it  Is  used. 


1.5  INSTALLATION  OF  NEW  TECHNOLOGY  MODULES 

Whenever  a  new  technology  module  Is  added  to  the  DYSCO  library,  certain  infor¬ 
mation  must  be  provided  to  the  executive.  This  Information  Is  provided  by 
"Installing"  the  technology  module  Into  the  common  block,  /XTM/,  which  con¬ 
tains  master  variables  and  arrays  describing  the  limitations  and  usages  of 
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that  technology  module.  In  addition,  FORTRAN  "calls"  are  inserted  into  speci¬ 
fic  predeveloped  subroutines  for  each  of  the  technical  modules  developed.  A 
new  shared  common  block  is  similarly  installed  via  the  common  block,  /XTM/, 
and  insertion  of  appropriate  "calls." 

1.6  THEORETICAL  BASIS 

The  basic  concepts  and  procedures  for  the  practical  use  of  DYSCO  have  been 
presented  in  this  general  discussion.  The  following  sections  include  detailed 
descriptions  of  component  and  force  modules  which  augment  material  presented 
in  Volume  II.  The  results  of  validation  analyses  for  advanced  helicopter  sim¬ 
ulations  are  presented  and  are  compared  with  results  from  other  sources. 

Except  where  noted  otherwise,  the  accuracy  of  the  analyses  is  only  dependent 
upon  the  limitations  of  the  approximations  made  by  the  user  in  the  development 
of  specific  component  representations  (data  sets). 

Volume  II  contains  a  complete  description  of  the  features  and  uses  of  the  pro¬ 
gram  and  the  necessary  information  for  the  development  and  installation  of  new 
components,  forces,  and  solutions.  Volume  III  contains  examples  of  input  and 
output  data  for  each  technology  module  and  several  modeling  examples. 


2.0  DYSCO  ANALYSES 


2.1  AOVANCED_ ELASTIC  BLADE  ANALYSIS  -  CRE3 

A  comprehensive  development  of  the  equations  of  motion  of  a  rotor  blade  was 
first  published  by  Houbolt  and  Brooks  (Reference  1)  in  1958.  These  equations 
were  reformulated  by  Hodges  and  Dowell  (Reference  2)  in  1974,  their  major 
contribution  being  improved  generality,  including  nonlinear  terms  and  the 
independent  verification  of  the  earlier  work.  Berman,  Giansante,  and 
Flannelly  (Reference  3)  expanded  these  previous  works  by  the  addition  of  five 
hub  degrees  of  freedom. 

The  present  research  uses  the  equations  of  motion  for  a  flexible  hingeless 
rotor  blade,  as  derived  in  Reference  2,  as  a  starting  point.  The  mathematical 
model  considered  for  the  rotor  blade  is  a  straight,  slender,  nonuniform  beam 
with  distributed  pretwist.  The  elastic  axis  (shear  center),  the  mass  center, 
and  the  area  centroid  are  considered  as  noncoincident.  Offset  at  r  =  e  is 
included  to  allow  for  consideration  of  the  hinge  offset  of  an  articulated 
rotor,  or  a  very  stiff  hub  of  a  hingeless  rotor.  Different  offsets,  ep  and 
e^,  are  used  for  flap  and  lead -lag  motion. 

The  rotor  blade  can  undergo  combined  In-plane  bending,  out-of-plane  bending, 
and  torsion,  and  the  rotor  hub  Is  assigned  six  degrees  of  freedom.  Hamilton's 
law  of  varying  action  (Reference  4)  has  been  used  to  derive  the  equations  of 
motion  in  a  generalized  coordinate  system  which  allows  a  direct  solution 
without  . the  consideration  of  force  equilibrium.  It  may  be  noted  that  the 
above  approach  is  equivalent  to  the  Lagrange  approach  if  the  shape  function  in 
time  domain  Is  not  being  used. 

2.1.1  Ordering  Scheme.  In  order  to  avoid  high  order  terms  that  might  overly 
complicate  the  equations  of  motion,  a  systematic  ordering  scheme  was  adopted 


to  determine  which  terms  should  be  Ignored.  The  ordering  scheme  employed  in 
this  study  for  the  various  parameters  Is  given  below: 


u/R 

-  0(£2) 

*?>  f/R  ■ 

0(e) 

v/R 

-  0(e) 

c,  t/R  - 

OU) 

w/R 

-  0(e) 

e/R  - 

OU) 

x/R 

-  0(1) 

^pc 

0(6) 

* 

-  0(e) 

XH’ 

V  ZH’  V 

<y  az  -  0(e) 

In  the  total  energy  equation,  if  any  terms  are  two  or  more  orders  higher  than 
the  lowest  order  term,  then  the  higher  order  terms  are  discarded. 

2.1.2  Coordinate  System  and  Motion  Variable.  As  demonstrated  in  Figure  1, 
the  coordinate  system  used  to  describe  the  blade  and  hub  motion  in  the  present 
analysis  is  as  follows: 

1.  Triplet  X,  Y,  Z  represents  an  inertia  frame,  R.  The  R  system 
coordinates  are  the  rotor  shaft  axes  when  there  is  no  hub 
motion. 

2.  Triplet  xp,  y,  z  is  fixed  in  a  reference  frame  Bj  which  rotates 
with  respect  to  R  frame  at  constant  angular  velocity  ft. 

3.  x,  y,  z  axes  are  also  fixed  in  B^,  in  which  x,  the  elastic  axis 
of  the  undeformed  beam,  is  inclined  to  xp  axis  at  the  precone 
angle  0pc.  This  system  is  referred  to  as  B  frame. 

4.  x',  y',  z'  are  a  blade  fixed  coordinate  system,  x'  is  tangent 
to  the  deformed  elastic  axis,  and  y',  z'  are  the  principal  axes 
of  local  cross  section.  Note  that  the  torsion  deformation  of 
the  blade  cross  section  is  about  the  x'  axis. 

The  deformations  of  elastic  axes  u,  v,  and  w  are  defined  with  respect  to  the  B 
system,  while  the  linear  and  angular  motion  of  the  rotor  hub,  X^,  Y^,  ZH,  <*x, 
oy  <*z  {Figure  2)  are  defined  with  respect  to  the  Inertia  frame,  R. 
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2.1.3  Hamilton's  Law  of  Varying  Action.  "Hamilton's  Principle,"  used  in  Ref¬ 
erence  2  to  derive  equations  of  motion  of  the  rotor  blade,  can  be  expressed 
as 


[5  (U  -  T)  -  5W]  dt  *  0 


(1) 


where  U  is  the  strain  energy,  T  is  the  kinetic  energy,  and  5W  is  the  virtual 
work  of  the  external  forces.  Hamilton's  law  of  varying  action,  which  is  a 
special  case  of  Hamilton's  Principle,  is  defined  by 


1 


[«  (U  T  T)  -  SW]  dt  -  2*-  Sq. 

3q,  ’ 


which  can  reduce  to  Equation  1  under  the  assumption 


£L 

aq, 


5q. 


i.e.,  "stationary"  in  calculus  of  variation's  terminology. 


The  Lagrange  equation 


(2) 


(3) 


d  ai_  3I_  ,  3U_ 

dt  3q. 


6W 

3q. 


5q.dt 


0 


(4) 


can  be  obtained  from  either  Equation  2  without  any  assumptions  or  Equation  1 
with  the  "stationary"  assumption. 


\ 
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It  should  be  emphasized  here  that  Hamilton's  law  of  varying  action,  Equation 
2,  is  not  merely  an  alternative  form  of  the  Lagrangian  equation,  Equation  4, 
but  is  much  more  general.  Equation  2  allows  one  to  pursue  the  solutions  of  a 
system  without  consideration  of  force  equilibrium  and  the  theory  of  partial 
differential  equations  (References  4  and  5),  while  Equation  4  does  not.  Un¬ 
fortunately,  a  detailed  discussion  of  the  differences  among  Equations  1,  2, 
and  4  is  beyond  the  scope  of  this  study. 

2.1.4  Strain  Energy  Contributions.  Using  the  derivations  of  Reference  2  with 
the  addition  of  root  springs,  the  strain  energy  contribution  can  be  expressed 
as 


5U 


{  Yv„6v"  +  Yy,6v'  +  Yw„6w"  +  Yw,6w'  +  Y0„60”  +  Y 


+  Y^  S<t>  +  Yu/6u'}  dx  + 


Kfw'6w'  +  K^'Sv'  +  KqW 


r  =  e 


(5) 


where  Kp,  KL,  and  are  flap,  lag,  and  torsion  spring  rate,  respectively,  and 


V  "  Mz'  cos^  +  *)  +  My'  sin(*  +  0) 

v  ■  V' 

Yw"  *  Mz'  s1n(*  +  t)  -  CO s(0  +  0) 

V  -  V*' 


V  “  Px' 

V  "  V  +  Tx' 

Y*  -  V"  -  Yw"v" 

V  -  V 


(6) 
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From  Equation  34  of  Reference  2 


V  /  a  O 

X'  XX 


,2  ,2 

d»7  df  «  EA  •  u'  +  ^  +  kA2  0'0'  -  eA  [v"  cos(0  +  0) 

* 


+  w"  sin(0  +  0)] 


V  -  L  ^xf  -  <V  d»  -  GJ*' 


V  ■  II  [9  +  t)'ovJn2  +  f2)  d*  df  -  EAk  ‘(0  +  #)'|u'  +  V  + 


v'  2.  w'  2 


2  T  2 


+  EB^'V  -  EB2*0'  (vn  cos  9  +  w"  sin  9) 


Aaxx  drj  df  »  ECj^"  +  ECj*  (w"  cos  $  -  v"  sin  $) 

J  J  A 


* 

My'  a  f*xx  d^  d^  *  EIy'tv"  s1n  (9  +  0)  -  w"  cos  (0  +  ^)]  -  ECj*0" 
A 


V  d?7  df  -  El  ,[v"  cos  [$  +  0)  +  w"  sin  (0  +  0)] 


f  V'2  w/2l 

-  EAeA  u  +  2~  +  -  EB2*0'0' 


\ 


From  Equation  35  of  Reference  2,  the  section  integrals  in  Equation  7  are 
defined  as  follows: 


A  *  JJA  dT?  df 

AeA  s  JJa  *» dr> 

« 

ly'  *  JJa  ^  d”  * 

V  s  jJA  *?2  d7?  df 

AIca2  *  JJa  (7?2  +  f2)  d,?  df 

J  a  JJa  (,?2  +  f2)  dr?  df 

Bl*  *  JJA  ( T)Z  +  r2)2  d7?  df 

B2*  s  JJA  +  f2)  dTj  df 

ci  s  JJa  *2  d7?  df 

ci*  s  JJa  d7?  df 

« 

2.1.5  Kinetic  Energy  Contribution.  The  position  vector  p  of  a  blade  element, 
after  deformation,  in  the  rotating  frame  B  is 

■♦444 

P  -  Xj  i  +  yj  j  +  Zj  k 


x  +  u  -  X0'  -  (v'  +  0$)  [t?  cos(0  +  0)  -  f  sin(0  +  0)] 


-  w'  [;?  sin(0  +  0)  +  f  cos(0  +  0) ]  r  ^ 


+  H 


+  H 


v  +  r0$  +  77  cos(0  +  0)  -  f  sin(0  +  0) 


w  +  r?  sin(0  +  0)  +  f  cos(0  +  0) 


(9) 
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The  rotational  velocity  of  frame  B  Is 


m 


where 


*z  +  Q 


cos  0pc  cos  #  cos  0pc  sin  #  sin  0pc 


-sin  } 


cos  i 


-sin  0pc  cos  i  -sin  0pc  sin  i>  cos  0pc 


(10a) 


The  translational  velocity  of  the  origin  of  frame  B  is 


-  m  TH 


Then,  the  total  velocity,  Vg  of  a  blade  element  can  be  expressed  as 

-+  +  -a  *  -a  -a- 

VB  *  VoB  +  *  +  WB  x  p  UZ) 

where  p  is  the  relative  velocity  of  a  blade  element  with  regard  to  rotating 
frame  B.  VQg  and  3g  x  p  are  contributions  from  the  motion  of  hub  degrees  of 


20 


freedom.  Thus,  from  Equations  9  and  11,  the  total  velocity  Vg  of  a  blade 
element  Is  given  by 


VB  "  V  +  Vyj  +  Vz  k 


(Xj  -  Oyj  cos  0pc  +  cos  0pc  cos  jfr  XH  +  cos  0pc  sin  0  YR 

•  4 

♦  s,n  V  ZH  +  “y  21  -  °z  >1*  1 

•  • 

+  (yj  +  QXj  cos  0pc  -  Ozj  sin  0pc  -  sin  $  XH  +  cos  j  YH 

-4 

-  nx  Zj  +  az  x j )  j 

•  • 

+  (Zj  +  flyj  sin  0pc  -  sin  0pc  cos  XR  -  sin  0pc  sin  i>  YH 


+  cos  etc  zh  +  nx  yj  -  fly  Xj)  k 


(13) 


where 


1 

n 

'  » 

•  j 

ft 

X 

X 

•°y 

'  -  [T]  • 

V 

n. 

•  ! 
a. 

z 

»  4 

z 

»  • 

(14) 


The  variation  of  kinetic  energy  is 
'R  rr 

- 


61  = 


Vg  •  6Vg  dm  dx 


(15) 
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Following  the  ordering  scheme  of  Reference  2,  with  the  additional  assumption 
that  the  order  of  hub  degrees  of  freedom  and  perturbation  shaft  speed  are 
0(e),  then,  from  Equations  13,  14,  15  and  the  application  of  Hamilton's  law, 
the  following  equation  can  be  obtained: 


rt. 


Tdt  -  6 q. 

1 


rR 


(Xu  «u  +  Xy  6m  +  Xw  6w  +  X0  64 


+  Xy,«  v'  +  Xw,  «w'  +  X^,  6<f>'  +  X^$  64$ 


+  \  SYH  +  XZH  «H  +  \  *x  +  X*Y  ‘*Y 


+  Xa  6 oz  Xx  5Xh)  dx  dt 
z  H 


(16) 


where 


[n2r  +  20  (v  +  r^$)j  -  m  sin  4  YH  -  m  cos  4  XH  +  mrflfl. 


Xu  -  m  irr  +  20 


(16a) 
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Xv  -  nfi2(v  +  r#$  +  e  cos  ($  +  *))  +  2mD(4pc  w  -  u) 

+  2me0(v'  cos  6  +  $  cos  6  +  w'  sin  6)  -  m(v  +  r0  ) 

* 

+  metf  sin  0  +  m(v  +  rtf  )  0  0  +  meO  0  cos  (0  +  0) 

d  tm 


K  16b) 


-  mQ  Ow  -  me  sin  6  0  0  +  m(v  +  rtf  )  0  0  +  me  cos  6  0  0 
y  y  s  z  z 

+  m  XH  sin  tf  -  m  Y^  cos  tf  +  mwflx  +  me  sin  6  0X+  mwOx 

+  mwfl^  -  mxO,  +  m  Zu  0  •  mr0R/.  0„  0  -  mrfl„ 
x  z  n  x  pc  y  x 

m/3pc  (O  r  +  20  (v  +  rtf$))  -  mw  -  me  cos  6  tf  -  mrflx  0 
■  m  Z||  •  si  (v  +  rtf$)  0X  +  mrfly 


•  • 

+  m  XH  sin  tf  (0)3pc  +  0X)  -  m  YR  cos  tf  (fl/3pc  +  flx) 

-  ""^z  Wpc  +  fix) 


r(  16c ) 


23 


X*  "  '  "V  *  '  m  qC  <  ^  '  Kmli>  cos  ^  +  sin  +  #> 


me02r  (*'  cos  *  -  v'  sin  6  -  sin  9) 


me  sin  9  Q2  (v  +  n>$)  -  meQ2  0pc  r  cos  0 


+  me  (v  sin  0  +  rrt  sin  6  -  w  cos  0) 


►(16d) 


2  2 

m^Km2  -  Kml  )  cos  6  sin  $  -  me  sin  6  sin  ^ 


+  me  sin  9  Y„  cos  f  +  me  sin  9  r  ft,  -  me  cos  9  r  0„  0 


-  me  cos  9  +  me  cos  9  r  Cl 


Xy/  »  -  me  (Q  r  cos  {9  +  4)  +  20  (v  +  )  cos  9) 


+  me  cos  9  XH  cos  1>  +  me  cos  9  Yu  sin  1>  -  me  cos  9  r  fin. 


Xw,  *  -  me  (0Z  r  sin  {9  +  *)  +  20  (v  +  )  sin  0) 


+  me  sin  9  XH  cos  ^  +  me  sin  9  Yu  sin  jfr  -  me  cos  9  r  00 


►(  16e ) 


-(16f) 


(leg) 


r  Xy  +  Xv, 


(16h) 
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Xo  "  dt  '  Vx  *sin  *  Z1  +  cos  *  sin  V  yl>  +  Vy  (cos  V  cos  *  Z1 


sin  ^pc  cos  *  ‘  vz  (cos  V  cos  ^  yl  +  sin  #  V  ' 


( 16  i ) 


X  ■  d- 
\  dt 


-Vx  (cos  \f>  Zj  +  sin  }  sin  B?c  yj)  +  Vy  (cos  0pc  sin  0  Zj 


+  sin  0pc  sin  ^  Xj)  +  Vz  (cos  0pc  sin  }  Xj  +  cos  0  Xj) 


(16j) 


x  *  d_ 
\  dt 


Vx  cos  *pc  *1  +  Vy  <sin  V  Z1  *  cos  V  V 


-  vz  sin  V  *if 


(16k) 


\  '  dt 


Vx  cos  V  cos  *  -  Vy  sin  i>  -  Vz  sin  Bpc  cos 


(161) 


YH  "dt 


Vx  cos  0pc  sin  +  Vy  cos  }  -  Vz  sin  Bpc  sin  f 


(16m) 


X  *  -  \ 

ZH  * 


V  sin  B  +  V  cos  B  r 
x  ^pc  z  Ppcr 


( 16n) 
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2.1.6  Virtual  Work  of  External  Forces.  The  virtual  work  for  generalized 
nonconservative  forces  may  be  expressed  as 


fR 


6W 


(Lu  5u  +  Ly  6v  +  Lw  6w  +  6j$)  dx  +  Lx^  6XH 


+  lYh  +  l2h  «H  +  \  Sax  +  \  4«y  +  \  iaz 


(17) 


where  L  and  M  are  generalized  forces  and  moments  for  each  degree  of  freedom. 

2.1.7  Governing  Equations  in  Physical  Coordinate  System.  Combining  Equations 
2,  5,  16,  and  17  yields  the  following  governing  equations: 

6u  Equation 


(V  «u'  -  Xu  Su  -  Lu  Su)  dx  -  0 


$y-Equati<?n 
fR 

(Yv.  5v"  +  Yy,  «v'  -  Xv,  5v'  -  Xv  8v  -  L,  6v)  dx 


(18) 


♦  Ikl  v'  ‘v'lp-e  '  ? 


(19) 


^w_Eqv.ati,on 

R 

(Yw"  5w"  +  V  5w'  ‘  Xw'  5w'  '  Xw  5w  '  Lw  5w)  dx 


+  [Kp  w'  «w']r.e  -  0 


(20) 
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Si  Equation 
'R 

(Y0"  +  V  W  +  Y*  *  •  X0  60  -  M0  60)  dx 


+  [K0  0  60]r-fi  -  0 


60  Equation 

X 


\  6*s  -  \  S*s 


dx  *  0 


6Xh  Equation 

fR 

Xy  dx  +  Lv  -  0 
0  aH  aH 


6Yh  Equation 
R 

XY  dx  +  Lv  *  0 
0  tH  yH 


6Zh  Equation 
R 

X7  dx  +  L7  -  0 
0  lW 


6a  Equation 
* 


X  dx  +  M. 
0  ax  0 


(21) 


(22) 


(23) 


(24) 


(25) 


(26) 
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&>v  Equation 


fR 


0  ay 


dx  +  M_  -  0 

y 


(27) 


Saz  Equation 


fR 


0  a2 


X^  dx  +  -  0 

az 


(28) 


If  hub  mass,  damping,  and  spring  rate  are  also  included,  then  Equations  23  - 
28  can  easily  be  rewritten  as 


SXH  Equation 


My  XH  +  CY  Xg  +  Ky  Xy 
n  A^j  n  n 


fR 

Xy  dX  +  Ly 

0  ah  ah 


(23a) 


6Yh  Equation 


V"  +  V«*V" 


Xy  dX  +  Ly 

0  th  th 


(24a) 


$ZH  Equation 


MZh  ZH  +  CZH  ZH  +  KZh  ZH 


fR 

X7  dx  +  L7 
0  £H 


(25a) 


6a  Equation 
—  ^ 


IA  a  +  C  ®Y  +  K  a 
\  x  ®x  x  \  x 


XA  dx  +  FT 
0  ax  °x 


(26a) 
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6a v  Equation 


X.  dx  +  M. 

0  ay  ay 


Saz  Equation 


'R 

0 


dx  +  Mft 


(27a) 


(28a) 


As  suggested  in  Reference  2,  the  tension,  T,  and  the  longitudinal  deflection, 
u,  can  be  eliminated  from  the  equations.  Considering  0  to  be  small,  with  0 
ignored  where  compared  to  0  in  the  nonlinear  terms,  the  equation  for  the  ten¬ 
sion  in  the  blade  becomes 


EA 


u' 


9'i' 


eA  (v"  cos  (0  +  0)  +  w"  sin  (0  +  0)H 


(29) 


Integrating  with  respect  to  x  and  solving  for  u  yields 


u'dx  - 


-  Ka2  0'0'  +  eA  (v"  cos  (0  +  0)  +  w"  sin  (0  +  0))  dx 


x  v'2  w/2 

<¥ +  ¥> dx 


(30) 


with  boundary  condition  u(0)  -  0. 
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Equations  31,  32,  and  an  expression  for  u'  developed  from  Equation  30  are 
substituted  into  Equations  19  -  28  and,  noting  that 

fR 

*  (  )  dx  ,  p  ■  me  sin  9  ,  and  q  -  me  cos  9  , 

J  Jo 

the  equations  of  motion  for  a  one-bladed  rotor  can  be  obtained  as 


30 


6v  Equation 


mv  Sv  - 


6v  + 


mx0s  Sv  - 


mXH  Sv  sin  $  + 


qXH  Sv'  cos  i> 


mY^  Sv  cos  0  - 


qYjj  Sv'  sin  0 


[p  +  mX/3pCj  <*x  Sv  cos  0 


(p  +  mX0pc]  ay  «v  sin  0  + 


mx<*z  Sv 


2qQv  Sv'  - 


2eA£JX  cos  0  Sv"  + 


2mfl  eA  v"  Sv  - 


2qfiv'  Sv 


2mfl  ( 


eA  sin  9  w")  Sv  - 


2pGw'  Sv 


2mfl  ( 


K,2  9'?)  Sv  + 


2qxfi  #s  Sv'  - 


2mQfl  w  Sv 
pc 


2qil  Sv 


2qxfi  az  6v'  - 


2qQ  az  6v 


T*v'  Sv'  + 


mQ  v  Sv  +  Kl  v'  Sv' 


r-e 


Eyw  w"  Sv"  + 


PfTtf  Sv 


pxfl  0  Sv'  + 


ECj*  sin  0  Sv"  + 


Ej  0'tf'  Sv" 


eA  T*  sin  6  <f>  Sv" 


+  (C0NST)v  +  (NONL IN) v 


Lv  6v 


Ev  v"  Sv" 


^2 


31 


(CONST) 


(NONLIN) 


33 


(Equation  34,  continued) 


+  -  eA  T*0  S w"  +  ECj*^"  Aw"  +  EjA'  9  f 


"  ' 

* 

meO  x  0  Aw' 


+  (C0NST)w  +  (N0NLIN)w  -  Lw  Aw 


where 


(CONST) 


w  -  -  eA  AT*  Aw"  +  mO*  A 


' 

__  x  Aw  +  TO2 
pc  , 


x  Aw' 


(NONLIN)  -  20 II  m  (v  +  i4  )  dx  w'  Aw'  -  e.  *  20 


ffR 

m  (v 

Jx 


+  r$$)  dx  6w" 


m  (v 


. 

m  (-Xh 

* 


+  m  (v  +  r^s)  Qx  Aw  +  m  (-X^  sin  )p  +  YH  cos  ^)  flx  Aw 


+  mx  a2  Qx  5w  +  m  (v  +  r$s)  QQy  Aw 


m  (XH  cos  i>  +  Yh  sin  ^)  0  Aw  +  m  (v  +  r$s)  ftx  Aw 


+  m  (v  +  r^$)  Qx  Aw  +  (EA  e£  A  sin  A  -  AE  cos  29)  v"  <f>  Av” 


+  w"  *  AE  sin  29  Aw* 


34 


84  Equation 


35 


where 


(CONST) 


K.2  9'  T*  64'  + 


*°2  l**  ■  S.12)  «"  • 


"°Vpc  X  cos  9  M 


(NONLIN) 


k.2  r 


20 


»  (v  +  r*$) 


-  EAK/  *' 


+  EAKa2  *'  eA  (v-  +  w "#) 


2eA° 


fR 


m  (v  +  r^s)  +  EA  eA  [ka2  9'  « 


-  eA  (v"  +  w"0)j 


(w"  -  v"#)  6$ 


AEI 


(w"2  -  v"2)  cos  9  sin  9  +  v"  w"  cos 


\ 


sin  9  69 


2 


29  \  69 
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£#  Equation 


37 


where 


(CONST) .  -  0 
"s 


(NONLIN), 


20 


mvdx 


<v'  +  60$) 


2m0 


[(v'  +  '#J  (v'  +  #c)  +  w'w']  dx 


r  St 


2m0  (v  +  n&$)  az  r6fs  + 


m  (XH  cos  i>  +  Yh  sin  p)  azr  60$ 


m  Ow  0y  rfiiSs  - 


mw  flx  r6fs 


2  mw  0X  r6fs  - 


mZH  fix  rS*s 


mx  0X  0y  r^$ 


38 


6Xh  Equation 


mv  sin  rf>  + 


qv'  cos  } 


pc 


w  cos  $  + 


pw'  cos  $ 


' 

•  < 

' 

p4  sin  $  - 

mx  1>s  sin  0  + 

4 

q  rf>$  cos  1> 

4 

+ 

m  +  m 

hJ 

P  + 


mx  p 


pc 


“y  • 


mx  az  sin  ^  + 


q  azcos  rp 


m  Qv  cos  $  +  2 


q  Qv'  sin  f 


pQw'  sin 


+  2 


mQ  /fpC  w  sin 


+  2 


pQ$  cos  0  + 


-  2 


mx  Q$$cos  ^ 


+  2 


qQ  $$  sin  0 


-  2 


s  y.  j 


mx  Q  az  cos  +  2 


qQ  az  sin  0 


2 

mQ  v  sin  rp  + 


2 

q  Q  v'  cos 


V w  + 


r  ^ 

pQ2w' 


COS  0 


* 

*  m 

* 

«« 

2 

mx  Q  rf>  sin  0  + 

J  4 

qQ2  cos  1> 

4 

+  ky 

H 

- 

(C0NST)y  + 
*H 


(N0NLIN)y  -  Ly 
*H  aH 


(37) 
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where 


J  (CONST)^ 
(NONLIN)y 


q  Irsln  $  + 


mxQ 


cos  + 
pc 


po2  « 


pc 


COS  i 


nix  (a*  cos  # 


f 


m  (v  +  n>s)  az  cos  ^ 


2 


*0  (v  +  n*s)  oz  sin  0  + 


mO  w  ax 


-  2 


mO  (v 


ris) 


az  Cos  if)  + 


m  w  a. 


y  2 


m  w  a. 


mx  oy  fly 
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6Yh  Equation 


* 

*  ' 

9 

9  ' 

mv  cos  }  - 

qv'  sin  $ 

* 

■Ope"* 

pw' 

sin  }  - 

- 

4 

V  j 

4 

* 

P0  cos  ^ 


9 

'  * 

9  *> 

* 

*  \ 

mx  cos  #  - 

J 

q  sin  j 

+ 

m  +  mv 

yhJ 

V 

« tfC  * 

• 

P 

• 

mx  az  cos  $  - 


2q  Qv'  cos  i> 


q  az  sin  f 


2p  flw'  + 


2m  flv  sin  $ 


2mO  0pc  ft 


cos  }  + 


2p  0$  sin  it 


2mx  flft$  sin  }  + 


2  qQaz  cos  ^ 


2  qfl  i>$  cos  0 


/•  r 


p 

mD  v  cos  f  + 


2  mx  fl  az  sin 


2 

q  0  v'  sin  0 


mftfl  w  + 
pc 


qfl  sin  j 


pfl2w' 


r  e 


sin  + 


mx  fl  cos  0 


+X  Yu 
/  Yh  h 


>  + 


( CONST )v  + 
YH 


(N0NLIN)v  -  Lv 
yH  yH 


41 


where 


(COHST)v 

yH 


(NONLIN)v  - 
TH 


q  Q^cos  #  + 


r  ^ 


•xO  cos  + 
pc 


PO2  B 


PC 


sin  ^ 


mx  (a^  sin  + 


«  (v  +  r*s)  az  sin  j 


-  2 


nfi  (v  +  n>$)  az  cos  + 


mO  w  a. 


-  2 


«fl  (v  +  r$$)  a2  sin  j  - 


m  w  a. 


-  2 


m  w  «x  + 


■K«x0y 
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6ZU  Equation 


raw  +  q  ^  +  m  +  Mz^  ZH  +  Mq  ax  cos  i>  +  Imx  ax  sin  } 


+  q  ay  sin  ^  -  mx  ay  cos  f  +  -  2q  Qa^sin  rf>  +  2mx  flax  cos 


*  2q«iy  cos  +  mx  0ay  sin  +  «z  ZH 


+  (CONST) 7  +  (N0NLIN)7  -  L7 

LW 


where 


(C0NST)7  -  0 


(NONL IN) 2  ■  m  (v  +  r^$)  fifly  +  2m  (v  +  r$$)  flx  +  2mx 
H  J 
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&av  Equation 


mx  4dc  v  cos  0  -  p v  cos  0  +  mxw  sin  0  +  qw  cos  0 
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Equation  41  (continued) 
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(NONLIN) 


x  cos  0  cos  0pc(NONLIN)y 


X  sin  ^  cos  0pC  (NONL IN) x 


in 


(v  +  n*s)  [yh  sin  +  +  XH  cos 


2m  (v  +  rtf  )  xfi  a, 
s  z 


2mO  (v  +  r^$)  (v  +  r$$) 


2.1.8  Ravleiqh-Ritz/Galerkin  Method.  According  to  the  Rayleigh-Ritz  method, 
arbitrary  functions  for  the  blade  elastic  displacement  can  be  separated  into  a 
sum  of  products  of  functions  of  r  and  t  only: 


v(r,t)  -  l  yi (t)  Yi (r)  »  I  Y1  Yi 
i 

(43) 

w(r,t)  -  l  Zj(t)  Zj(r)  -  l  Zj  Zj 

(44) 

j 

0(t , r )  -  l  *k(t)  ♦k(r )  -  l  tfk  #k 

(45) 

k 


where  Y..(r),  Zj(r),  *k(r)  are  modal  functions  and  y^,  Zy  tfk  are  generalized 
coordinates.  Substituting  these  sums  into  Equations  33  -  42  yields  the  equa¬ 
tions  of  motion  for  all  the  generalized  degrees  of  freedom.  These  can  be 
automatically  handled  by  the  computer  program.  Therefore,  Equations  33  -  42 
represent  the  final  form  of  the  equations  of  motion  and  no  further  derivation 
is  necessary.  Note  that  the  modal  functions  are  chosen  as  any  complete  set  of 
linear  independent  functions  that  satisfy  all  the  geometry  boundary  conditions 
(but  not  necessarily  the  force  boundary  conditions),  and  are  piecewise  twice 
differentiable. 


It  can  be  shown  that  this  Rayleigh  Ritz  approach  is  equivalent  to  the  Galerkin 
method  as  applied  to  equations  of  motion,  including  all  the  "trailing  terms" 
(Reference  6).  Note,  however,  that  If  the  Galerkin  method  applies  to  equa¬ 
tions  of  motion  only  and  does  not  include  the  "trailing  terms,"  then  the  modal 
functions  have  to  satisfy  all  geometry  and  force  boundary  conditions.  This  is 
not  always  possible  when  masses  or  dampers  are  at  the  boundaries.  In  addition 
to  this,  the  modal  functions  have  to  be  differentiable  to  the  degree  required 
by  the  differential  equations  (for  the  present  study,  fourth  differentiable), 
which  is  two  degrees  higher  than  the  usage  of  energy  method. 

2.1.9  Program  Features  of  CRE3.  Three  physical  degrees  of  freedom  are 
assigned  to  each  blade:  out-of-plane  bending,  In-plane  bending,  and  torsion. 
The  rotor  speed  perturbation  Is  treated  as  an  additional  rotor  degree  of  free¬ 
dom  which  may  couple  with  transmission  and  fuel  systems  to  reflect  the  feed¬ 
back  relationship  between  power  supply  and  rotor  speed.  However,  component 
representations  of  these  systems  are  not  yet  available.  All  four  blade  physi¬ 
cal  degrees  of  freedom  mentioned  above  are  optional  in  CRE3,  except  that  at 
least  one  blade  elastic  degree  of  freedom  must  be  chosen.  The  bending  and 
torsion  motion  of  the  blade  Is  expanded  as  series  of  the  product  of  decoupled 
modes  and  generalized  degrees  of  freedom  (as  in  Equations  43,  44,  and  45). 

Up  to  five  modes  for  out-of-plane  bending  and  three  modes  for  in-plane  bending 
and  torsion  may  be  used.  The  user  has  the  option  of  either  inputting  the  mode 
shape  and  its  first  and  second  derivatives,  or  letting  the  program  generate 
the  mode  shape  automatically  based  on  the  normal  modes  of  a  nonrotating  beam 
and  the  user-selected  boundary  conditions.  Because  the  derivation  of  equa¬ 
tions  of  motion  in  this  present  study  is  based  on  the  energy  method,  the  sat¬ 
isfaction  of  geometric  boundary  conditions  for  each  mode  shape  may  not  be  re¬ 
quired,  although  it  is  always  a  plus  if  all  boundary  conditions  are  satisfied. 

Theoretically,  a  mode  shape  may  be  obtained  by  Integrating  its  2nd  derivative 
twice  while  applying  the  appropriate  displacement  and  slope  boundary  condi¬ 
tions.  However,  for  a  discretized  system*  this  approach  results  in  numerical 
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problems,  especially  for  higher  order  modes.  So,  a  more  complicated  input 
procedure  has  been  adopted.  The  mode  shape  and  its  first  and  second  deriva¬ 
tives  must  be  input,  instead  of  the  second  derivative  and  its  associated  bound¬ 
ary  condition. 

The  mode  shape  data  are  defined  at  blade  stations,  starting  at  the  blade  root 
and  ending  at  the  tip.  Hinge  and  bearing  locations  are  input  as  indices  of 
the  blade  station.  Based  on  the  derivation  of  Reference  2,  the  blade  rotation 
order  is  in-plane,  out-of-plane,  and  torsion.  Care  must  be  taken  when  simu¬ 
lating  an  articulated  rotor  to  ensure  that  hinge  and  bearing  locations  are 
consistent  with  the  blade  rotation  order. 

In  CRE3,  six  hub  degrees  of  freedom  in  the  fixed  system  (none  of  which  is  man¬ 
datory)  are  assigned  to  the  hub.  In  a  model,  these  hub  degrees  of  freedom  may 
automatically  couple  with  a  fuselage  component  (CFM2)  if  the  component  number 
of  the  rotor  coincides  with  any  one  of  the  assigned  rotor  numbers  in  component 
CFM2.  The  hub  degrees  of  freedom  may  also  couple  with  any  degree  of  freedom 
with  the  same  name  through  the  OYSCO  naming  convention. 

CRE3  can  couple  with  the  control  system  (CCE0,  CCE1)  only  if  the  torsional 
degree  of  freedom  is  activated.  In  that  case,  pitch  horn  length  and  pitch 
horn  station  data  will  be  requested.  If  a  model  is  formed  by  including  both 
the  rotor  component  and  control  system,  the  same  component  number  has  to  be 
assigned  to  each  in  order  to  successfully  couple  them.  In  addition  to  coup¬ 
ling  with  the  control  system,  the  blade  may  also  couple  with  CSF1,  CFM2,  or 
CES1  via  the  OYSCO  naming  convention,  if  the  blades'  implicit  degrees  of  free¬ 
dom  are  chosen,  or  through  the  use  of  component  CLC1.  Although  programmati¬ 
cally  OYSCO  will  perform  the  coupling,  logically  the  mass  coupling  of  these 
components  with  the  rotor  will  not  be  correct  unless  the  component  mass  mat¬ 
rices  are  derived  for  a  rotating  system. 

In  CRE3,  if  a  term  in  the  mass  matrix  depends  on  the  acceleration  of  a  degree 
of  freedom,  the  delay  acceleration  of  the  DOF  is  used  in  order  to  avoid  any 

iterative  procedure  in  the  solution  method.  This  approach  is  justifiable  if 

v 

■j 
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the  nonlinear  term  is  not  the  lowest  order  term  in  an  equation.  This  is  the 
case  in  the  present  study. 

CRE3  may  also  couple  with  any  one  of  the  aerodynamic  force  modules  FRA0,  FRA2 
and  FRA3.  FRA0  is  the  simplest  aerodynamic  force  module  and  deals  only  with 
the  linear  aerodynamics.  FRA2  can  handle  airfoil  tables,  and  can  read  the 
inflow  distribution  from  an  outside  sequential  file.  Both  FRA0  and  FRA2  con¬ 
sider  only  steady  aerodynamics  without  yaw  flow  correction.  The  C81  type  of 
aerodynamics  is  installed  in  FRA3.  A  detailed  description  of  FRA3  is  given  in 
the  next  section.  In  addition  to  aerodynamic  force  modules,  CRE3  may  also 
couple  with  the  vibration  shaker,  FSS1,  to  get  useful  information  about  the 
rotor  under  harmonic  excitation. 

2.2  ADVANCED  ROTOR  AERODYNAMIC  ANALYSIS  -  FRA3 

The  basic  rotor  aerodynamic  analysis  used  in  the  present  study  is  based  on  the 
report  by  McLarty  (Reference  7).  The  steady  state  lifting  coefficient  CL, 
drag  coefficient  Cq,  and  moment  coefficient  C^  can  be  either  calculated  or 
looked  up  in  a  table.  For  either  option,  the  angle  of  attack  and  Mach  number, 
based  on  the  two-dimensional  strip  theory,  are  modified  by  the  yaw  flow  angle 
to  account  for  the  three-dimensional  aerodynamic  effect.  The  analysis  and 
experimental  work  presented  in  References  8  and  9  are  used  to  find  the  induced 
velocity  distribution  across  the  rotor  disk. 

The  unsteady  aerodynamic  effects  on  the  pitching  moment  coefficient  (Refer¬ 
ence  10),  on  the  lift  coefficients  (References  11  and  12),  and  on  the  drag 
coefficient  (Reference  12)  are  included  to  modify  the  steady  lift,  drag,  and 
moment  coefficients.  Figure  3  shows  the  flowchart  for  the  rotor  aerodynamics 
logic.  The  equations  in  the  following  sections  can  be  found  in  Reference  7 
but  are  reproduced  in  this  report  for  the  sake  of  clarity. 

\ 

2.2.1  Induced  Velocity  by  Eauation.  The  theoretical  basis  of  the  induced 
velocity  calculation  in  this  paragraph  is  given  in  Reference  9.  Several 
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Figure  3.  Rotor  Aerodynamic  Logic. 
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assumptions  of  the  vortex  theory  adopted  In  Reference  9  are: 

1.  The  contraction  of  the  slipstream  Is  neglected. 

2.  The  number  of  rotor  blades  Is  Infinite  and  the  slipstream  does 
not  rotate. 

3.  The  circulation  Is  constant  along  the  blade  radius. 

4.  The  angle  between  the  slipstream  and  the  rotor  disc  at  each 
point  of  the  rotor  disc  is  constant. 

Under  the  above  assumptions,  the  induced  velocity  across  the  rotor  is  then 
computed  from  the  contribution  of  bound,  free  tip,  and  shed  vortices.  Like 
the  simple  momentum  theory,  this  approach  fails  when  n  (advance  ratio)  and  X 
(inflow  ratio)  are  both  equal  to  zero.  In  that  case,  the  concept  of  a  solid 
disk  is  used  to  replace  a  rotor  to  make  the  result  realistic.  Further  refine¬ 
ment  of  this  theory  (Reference  8)  permits  its  application  to  a  wide  range  of 
rotor  operating  states  such  as  autorotation,  vortex  ring,  windmill  state, 
climbing,  and  high  speed  flight. 

Assuming  that  the  rotor- induced  velocity  is  parallel  to  the  centerline  of  the 
rotor  shaft,  the  local  induced  velocity,  v ,  can  be  expressed  as 

v  -  v  Fn  (46) 


where 


v  ■  the  average  induced  velocity  across  the  rotor  disk 
Fjj  -  an  induced  velocity  distribution  function. 

Fn  is  a  function  of  the  advance  ratio,  p,  radial  station,  x,  blade  azimuth, 
and  inflow  ratio,  X.  It  may  be  either  calculated  or  looked  up  in  a  data  table 
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(see  paragraph  2.2.5).  The  equation  used  to  represent  is 


FN  -  f  x  [i  +  fj(/i)  cos  tf] 


+  f2(x,^)  fj(/i)  K27  /-.5  V*  +  /25  vj  +  (*/•)£ 


where 


i/  is  the  average  induced  velocity  across  the  rotor  disc 
x  Is  the  nondimensional  station 
$  Is  the  blade  azimuth  angle 
K27  Is  the  tip  vortex  coefficient 

VN  Is  the  flight  path  airspeed  in  ft/s  divided  by  1.0  ft/s 
(v)N  is  v  in  ft/s  divided  by  1.0  ft/s 
and  f j(/x)  is  defined  as 


Mm) 


.5,  if  0  <  1  rad/s 

11.25  n,  if  0  >  1  and  m  <  -1067 

1.36  -  1.5  m,  if  Q  >  1  and  .1067  <  v  <  .5733 
.5,  if  Q  >  1  and  n  >  .573 


Function  f2  (x,tf)  is  given  by 


f2(x,0) 


0,  if  (x  >  .7) 

or  (105*  <  *  <  255*) 
or  (315*  <  0  <  360*) 
or  (0*  <  j  <  45*) 

sin  [6(tf  -  45*)],  If  (x  >  .7)  and 

if  (45*  <  x  <  105*) 
or  [(225*  <  i>  <  351*)] 


The  average  Induced  veloqity/is  given  as 

, _ ?B_ 


L+8x|*l/3) 

(CB  +  8  \£)  (CB  +  8  mZ) 


(47) 


(48) 


(49) 


(50) 
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Cg  Is  defined  as 

CB  •  tyin2  -  xj) 

where 

Cy  is  the  thrust  coefficient 
Xh  is  the  hub  extent  divided  by  rotor  radius 
B  is  tip  loss  factor. 


/ 

/  / 


The  ground  effect  is  included  if  the  center  of  gravity  of  the  vehicle  is 
between  one-fourth  and  one  rotor  diameter  above  the  ground  and  the  airspeed 
is  less  than  30  ft/s: 

v*  -  v  [1  +  (G  -  1)  ((v  -  30)/30)2]  (52) 

where 

_* 

v  is  the  induced  velocity  with  the  ground  effect  included 
G  is  0.25,  plus  the  altitude  of  the  CG  divided  by  the  rotor  diameter 
V  is  the  flight  path  velocity  in  ft/s. 


2.2.2  Steady  State  Aerodynamics.  If  the  relative  velocity  of  a  blade  element 
with  regard  to  the  air  Is  the  vector  sum  of  tangential  velocity  Uy  (in  the  hub 
plane,  positive  in  the  blade  drag  direction),  the  vertical  velocity  Up  (per¬ 
pendicular  to  the  hub  plane,  positive  downward),  and  the  radial  velocity  l)R 
(positive  outward),  then  the  standard  two-dimensional  strip  theory  gives  us 
the  inflow  angle  4^  defined  by 

-1 

*A  *  tan  ^  (53) 

and  the  angle  of  attack 

a  =  *A  +  4  +  0A  (54) 

where 

4  is  the  torsion  deformation  and  $ ^  is  the  sum  of  the  root 
pitch  angle  and  built-in  twist. 
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The  three-dimensional  effect  of  aerodynamics  may  be  included  by  the  consider¬ 
ation  of  yaw  flow  angle,  A,  which  is  defined  as 


A 


tan 


-1 


He 


(55) 


If  |a|  >  60*,  then  A  is  given  by 


A  ■  60*  •  sin  (A)  (56) 

Then  the  modified  angle  of  attack,  arm0(j,  and  the  effective  Mach  number,  M^, 
are  expressed  as 


a 

M 


mod 

eff 


a  cos  A 


UR  +  UT  +  UP 
u 

sound 


(cos  qj  A) 


«2 


(57) 

(58) 


As  suggested  by  References  12  and  13,  Equations  55  -  58  should  be  adopted  in 
all  aspects  of  the  steady  state  lift  determination.  Both  qj  and  q2  in  Equa¬ 
tion  58  must  be  supplied  by  the  user  and  should  be  determined  carefully.  Note 
that  Hoerner's  result  (Reference  13)  can  be  approximated  by  assuming  qj  =  2, 
Q2  m  1  dj  “  1*  ^2  "  *5* 


The  steady  state  lift  coefficient  is  then  expressed  as 


CL$  *  CL  (fltmod’  Meff)/cos  A 

Equation  59  is  valid  whether  is  derived  or  looked  up  in  a  table. 


(59) 


According  to  Harris  (Reference  12),  the  blade  drag  coefficient  in  the  blade 
element  plane  is 

CDn  •  CD  Meff>  (60) 
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and  the  drag  coefficient  in  the  radial  direction,  based  on  U^,  is  given  by 

CD  -  CD  (a  -  0,  M  -  .3)  (61) 

R 

Equation  61  is  used  in  addition  to  Equation  60  because  the  resultant  drag 
force  should  be  in  the  same  direction  as  the  total  wind  velocity. 

The  steady  state  moment  coefficient  is 

CHS  ■  CH  <«•  "eff>  <«> 

The  program  provides  the  option  of  using  either  formulas  or  lookup  tables  to 
determine  the  steady  state  coefficients  C^,  Cp,  and  C^.  Details  of  the  for¬ 
mulas  are  given  in  paragraph  2.2.4. 

2.2.3  BUNS  Unsteady  Aerodynamic  Model.  As  stated  in  References  10,  11,  and 
12,  the  assumption  of  quasi -steady  aerodynamics,  which  ignores  the  time  rate 
of  change  of  the  blade  angle  of  attack  and  the  virtual  mass  effect  of  the  air, 
is  too  crude  to  be  adequate  in  certain  situations.  Also,  it  is  impractical  to 
obtain  meaningful  wind  tunnel  test  data  over  the  full  range  of  angles  of 
attack  and  yaw  angles,  or  to  construct  a  theoretically  general  unsteady  aero¬ 
dynamic  model.  Therefore,  a  procedure  based  on  limited  wind  tunnel  data  and  a 
simplified,  unsteady  two-dimensional  theory  is  used  to  determine  the  unsteady 
aerodynamic  coefficients  as  described  in  the  following  paragraphs. 

2.2.3. 1  Pitching  Moment  Coefficient  -  The  determination  of  pitching  moment 
coefficient  is  based  on  the  work  of  Carta,  et  al .  (Reference  10).  (The  theo¬ 
retical  basis  for  their  unsteady  aerodynamics  can  be  found  in  Reference  14.) 
Their  test  data  was  obtained  from  differential  pressure  transducers  mounted  on 
a  2-ft-chord  NACA  0012  model.  Testing  was  performed  in  the  two-dimensional 
channel  of  the  UARl  8-ft  wind  tunnel. 

Carta,  et  al.,  assume  that  "...sinusoidal  data  could  be  generalized,  through 
cross  plots,  to  functions  of  instantaneous  angle  of  attack,  angular  velocity 
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parameter.  A,  and  angular  acceleration  parameter,  B,  for  a  given  Mach  number." 
From  Reference  10,  A  and  B  are  given  by 

A  -  (fy)  a  (63) 

8  ■  t|y)2  i  (64) 

where  c  is  chord  length',  U'is  wind  velocity  perpendicular  to  airfoil  leading 
edge,  a  is  time  rate  change  of  angle  of  attack,  and  <f>  is  the  blade  torsion  ac¬ 
celeration.  Carta's  table,  adopted  here,  excludes  the  contribution  from  the 
steady  state  aerodynamics,  i.e.,  ACM  -  0  when  A  »  B  -  0,  and  is  based  on  Mach 
number  equal  to  zero.  The,  effect  of  Mach  number  on  the  stall  point  is  in¬ 
cluded  by  computing  a  shift  in  the  angle  of  attack  argument  before  entering 


the  A  -  B  tables  as  follows  (Reference  6): 

“Carta  "  1*1  for  Meff  *  2  (65) 

“Carta  *  1*1  13.5  -  16!is5(Meff  -  .2)  for  2  -  Meff  <  *6  (66) 

*Carta  *  1*1  *  1,93  for  Meff  -  ,6  <67> 

Then  AC^  is  determined  as 

ACm  -  (sin  a)  AC„  (A,  B,  <»Carta)/N  (68) 

where  N  is  used  to  consider  the  compressibility  effect: 

N  *  1 1  *  Meff2  for  Meff  <  '6  (69) 

N  -  .8  for  Meff  >  .6  (70) 
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The  total  pitching  moment  coefficient,  C^,  can  then  be  expressed  as 

CM  "  CM  <steady  state)  +  ACm  (71) 

2. 2.3. 2  lift  Coefficient  -  The  lift  coefficient  determination  is  based  on  the 
work  of  Scanlan  and  Rosenbaum  (Reference  11)  and  the  assumption  that  the  fre¬ 
quency  of  oscillation  of  the  rotor  blade  is  small  (i.e.,  the  value  of  circula¬ 
tion  function  is  equal  to  1).  The  lift  coefficient  increment  of  a  flat  plate 
airfoil  due  to  the  effect  of  unsteady  aerodynamics  in  an  unstalled  region  is 


AC, 


2* 


.  fci-  _  /a  _  lx  i£  .  btt 

,„2  o,i2  la  Z]  U  +  2U 


2U  2U 


(72) 


where 


z  is  vertical  acceleration 
d  is  torsion  acceleration 
$  is  torsion  velocity 
U  is  the  resultant  wind  velocity 
a  is  angle  of  attack  change  rate 
b  is  the  semi chord 

a  is  the  distance  from  midchord  to  the  elastic  axis,  divided  by  b. 


The  first  and  second  terms  on  the  right  side  of  Equation  72  represent  the  vir¬ 
tual  mass  effect  of  the  air  due  to  the  acceleration  of  the  rotor  blade  in  the 
vertical  and  torsion  directions,  respectively.  The  third  term  is  the  effect 
of  blade  torsion  velocity  on  the  angle  of  attack.  The  last  term  is  a  damping 
term  based  on  rate  of  change  of  the  angle  of  attack  a.  If  the  elastic  axis  is 
located  at  the  quarter  chord,  the  value  of  a  in  Equation  72  is  minus  one-half. 


As  suggested  in  Reference  12,  the  dynamic  stall  effect  should  be  included  as 
follows: 

For  “mod  >  “stall 
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V  *  "mod  •  (s 


* 

ign  a)  Min  « 

» 


where 


La  -  61.5  In  (.6/MM)  J\a  b/U| 


.3  for  Mgff  <  .3 
Meff  for  .3  <  Meff  <  .6 
.6  for  Me^  >  .6 


If  “mod  -  “stall 


ttRL  "  “mod 


The  total  lifting  coefficient  is  given  by 


CU  *  CL  (aRL’  Meff*  aD,  +  ACL 


where  CL  (aR^,  is  the  lift  coefficient  corresponding  to  arRL  and  M0^ 

under  stationary  conditions. 


2. 2.3.3 
fied  as 


-  According  to  Reference  12,  the  Cn  should  be  modi 

UN 


CDn  *  CD  (ttRD’  Meff> 


where 

“RO  '  «  -  <si"  •>  Min  '  f  z  «  |  ' 

.  4 


2.2.4  Steady  State  Aerodynamics  bv  Equations.  It  is  noted  that  all  the  equa¬ 
tions  in  this  paragraph  are  reproduced  from  Reference  6  for  the  sake  of  com¬ 
pleteness.  In  order  to  facilitate  the  derivation  which  follows,  an  array  Y  is 
first  defined  as  shown  below. 


61 


Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 


1 


(7) 

(8) 

(9)  } 

(10U 

(11) J 

(12) 

(13) 1 

(14) J 

(15) 

(16) 

(17) 

(18) 

(19) 

(20) 
(21) 


Drag  divergence  Mach  number  for  a  -  0 

Mach  number  for  lower  boundary  of  supersonic  region 

Maximum  C^,  normal  flow,  M  *  0 

Coefficients  of  Mach  number  in  maximum  C.  equation, 
normal  flow  L 


Maximum  CL,  reversed  flow,  M  =  0 
Slope  of  lift  curve  for  M  -  0 

Coefficients  of  Mach  number  for  lift  curve  slope 
subsonic  region 

CQ  for  a  -  0,  M  *  0 

Coefficients  of  a  in  nondivergent  drag  equation 


(/deg) 


in 


(/deg) 

(/deg) 

(/deg) 


(/deg) 


(/deg) 

(/deg2) 


Coefficient  In  supersonic  drag  equation 
Maximum  nondivergent  CQ 
Thickness/chord  ratio 


Not  used 


Drag  rise  coefficient  (/deg) 

Coefficient  of  yaw  angle  in  Mach  number  equation 
Exponent  in  Mach  number  equation  for  yawed  flow 


m  2 

Y  ^22'  (  Coefficients  of  a  for  Mach  critical  in  steady  CM  (/de9  ) 

Y  (23)  equation  "  (/deg) 

Y  (24) J 

Y  (25)  CM  for  a  -  0,  M  -  0 


2.2.4. 1  Lift  and  Drag  Coefficient  -  The  slopes  of  lift  curve  in  subsonic, 
transonic,  and  supersonic  flow  regions  are  defined  as 

*  Y(8)  +  Y( 9)  M  +  Y ( 10)  M2  +  Y(ll)  M3  (subsonic)  (79) 

a2  -  Bq  +  Bi  M  +  B2  M2 


(transonic) 


(80) 


a3  -  4/(57 .3  MZ  -  1) 


(supersonic) 


(81) 


where  M  is  the  effective  Mach  number  defined  by  Equation  58.  The  lift  curve 
slope  corresponding  to  the  critical  Mach  number  (in  the  subsonic  region)  is 


-  a,  (Y(l)) 


(82) 


The  Mach  number,  MjC,  used  in  calculating  the  lift  curve  slope  (in  the  super¬ 
sonic  region)  can  be  found  from  Equation  81  as 


"sc  -  A  +  ((• 


0698)/(a,)  ) 


CR 


(83) 


Then,  the  Mach  number  of  the  transonic  region  is  defined  between  Y(l)  and 
where 


M$  -  Max  (Y (2) ,  M$c} 


(84) 


The  three  coefficients  Bg,  Bp  Bg  in  the  transonic  region  are  determined  by 
the  following  three  equatipns: 


(V(l))  -  a,  (V(l))  / 
a2  <MS>  *  a3  <MS> 

d&A 

dH  <MS>  ’  dM  <MS> 

/  / 

The  modified  angle  of  attafck  is  given  as 
or  cos  A  if  | or |  <  90* 

a  *  * 

1  [a  if  | a |  >  90’ 

The  typical  curve  of  versus  or  is  shown  in  Figure  4. 


(85) 

(86) 

(87) 


(88) 
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Figure  4.  General  Lift  Coefficient  Versus  Angle  of  Attack  Curve. 
In  the  normal  flow  region,  |a-j|  <  90_: 
or  «  |<*j| 

CL  -  Y(3)  +  Y(4)  M  +  Y(5)  MZ  +  Y(6)  M3 
ar$  *  C^/a 

aB  *  ®S  +  5- 

In  the  reverse  flow  region,  |a|  >  90*: 

a  -  180*  -  |a1 | 

CL  -  Y(7) 

®B  *  CL/a  +5* 

For  «B  <  a  <  90* 

CL  »  [{1.876  sin  a  -  (.581))  K  +  .81]  cos  a 
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where 


\ 


’1  +  .25  M4  ,  if  H  <  1 

K  »  ■ 

..85  +  .82/(M  -  .8),  If  H  >  1 

The  curve  of  Cg  against  a  is  seen  in  Figure  5.  The  corresponding  Cg  and  a  at 

P3  are  CDX  and  °X’  resPectively*  At  P3,  either  ax  *  a$  and  CQX  <  Y(16),  or 
ax  <  as  and  CQX  ■  Y( 16) . 


Figure  5.  General  Drag  Coefficient  Versus  Angle  of  Attack  Curve. 

Unlike  the  lift  coefficient,  determination  of  the  drag  coefficient  requires 
the  definition  of  drag  characteristics  in  only  two  regions:  the  subsonic 
and  the  supersonic.  In  the  subsonic  region,  CQ  is  expressed  as 

For  ax  >  o  >  0 

CQ  -  Min  { Y( 16) ,  (Y(12)  +  Y(13)a  +  Y(14)a2  +  Max  {0,  Y(19)a 

-  Y(l)  +  Max  (M,  .35)))) 
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For  90*  >  a  >  a* 

2  2 

CD  *  *4  s*n  a  *  (CDX  '  *4  s*n  aX^  cos  ®/cos  ax 

where 

K4  -  2.1  K 

In  the  supersonic  region,  M  >  Mg,  Cq  Is  the  function  of  both  a  and  M  and  Is 
given  as 


C 


0 


Min 


Y ( 16) ,  [y(12)  +  4  [(a/57.3)2  +  Y(15)]/v4l2  -  l]  - 


2. 2. 4. 2  Moment  Coefficient  -  The  moment  coefficient,  CM,  Is  a  function  of 
angle  of  attack,  a,  effective  Mach  number,  Mg^,  and  some  Input  constants 
which  define  the  characteristics  of  the  CM  curve.  The  Mach  number  of  the 
break  point  can  be  expressed  as 

"b  ■  *3  ♦  A2  M  ♦  A1  l«|2 

-  Y(24)  *  Y(23)  |o|  +  Y(22)  |«|2 

which  Is  defined  as  Mach  number  at  the  point  where  CM  breaks  rapidly  away  from 
an  input  constant  or  Y(25).  The  angle  of  attack  corresponding  to  a  partic¬ 
ular  Mb  In  this  equation  Is  ag.  The  value  of  ag  at  Mb  -  0  is  defined  as  A&. 
Then,  the  determination  of  Is  described  as  below: 

For  ag  >  a 

CM  -  Y(25) 
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For  A5  >  o  >  aB 

CM  -  (|a|  -  oB)  3  sign  (a)  +  Y(25) 

where  3  depends  on  and  an  input  critical  value,  Y(24),  which  is  the  point 
on  the  break  curve  for  a  =  0: 


.045  +  .0425  Meff 

for  Meff  *  A3 

.0!  -  .2  <Heff  -  Aj) 

for  Meff  >  A3 

For  |ct|  >  A5 

CH  *  (A5  ‘  3  ‘  '°0646  < lal  *  ty  si9n  (®)  +  Y (25) 

For  | a |  >  90* 

CM  »  -  .5  (CL  cos  a  +  Cp  sin  a)  +  Y(25) 

It  is  assumed  that  the  aerodynamic  center  is  located  at  the  .75  chord,  rather 
than  at  the  .25  chord.  Under  this  assumption,  the  pitching  moment  is  mainly 
due  to  lift  and  drag  forces. 

2.2.5  Induced  Velocity  bv  Data  Table.  When  a  data  table  is  used,  the  local 
induced  velocity,  »/.,  can  be  expressed  as  a  function  of  advance  ratio,  /z, 
average  inflow  ratio,  radial  blade  station,  r/R,  and  blade  azimuth  angle, 


^  (M,A,r/R,0) 


NH 

a0  +  l  an  cos  (ntf)  +  bn  sin  (ntf) 
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where 


Is  the  average  induced  velocity  across  the  rotor  disc 

a  .  are  coefficients  of  the  harmonics 
n  n 

NH  is  the  highest  harmonic. 

For  a  data  table,  the  permissible  number  of  advance  ratios  is  3,  the  number  of 
inflow  ratios  is  2,  and  the  number  of  harmonics  is  6.  If  n,  or  A,  or  r/R  is 
not  the  same  as  used  in  FRA3A,  then  subroutine  WKFIND  interpolates  the  values 
wherever  the  computed  values  are  within  the  range  of  their  respective  input. 
If  the  computed  value  is  outside  the  range  of  its  input,  the  procedure  uses 
the  input  value  which  is  closest  to  the  computed  value. 

2.2.6  Program  Features  of  FRA3.  In  FRA3,  the  user  has  the  option  to  choose 
whether  the  steady  state  aerodynamic  coefficients  are  to  be  calculated  by 
equation  or  through  table  look-up.  The  basic  independent  variables  used  by 
both  the  equations  and  the  table  look-up  procedure  are  angle  of  attack  and 
Mach  number.  As  shown  in  the  preceding  section,  the  angle  of  attack  and  Mach 
number  are  modified  based  on  yaw  flow  angle,  furthermore,  the  blade  drag  in 
the  radial  direction  Is  also  considered. 

The  BUN's  unsteady  aerodynamics  described  by  Equations  63  -  78  are  also  imple¬ 
mented  in  FRA3  as  an  option  under  the  assumption  that  the  rotor  reduced  fre¬ 
quency  Is  small.  It  is  emphasized  here  to  caution  the  user  that,  according  to 
Reference  7,  "...for  rotors  with  cambered  airfoils  where  the  chord  line  and 
zero  lift  line  are  not  coincident,  it  is  advisable  to  use  data  table  rather 
than  equations  to  compute  the  aerodynamic  coefficients,"  because  the  original 
equations  were  developed  based  on  the  symmetric  airfoil.  Note  that  both 
options  take  the  effect  of  the  yaw  flow  into  consideration  and  include  the 
blade  drag  In  the  radial  direction. 

The  same  options  (equation  or  data  table)  are  also  available  for  the  rotor 
Induced  velocity  calculation.  According  to  Reference  7,  the  induced  velocity 
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equations  described  In  Equations  46  -  52  are  supposed  to  be  valid  for  a  very 
wide  range  of  rotor  operating  states,  including  autorotation,  vortex  ring, 
windmill  state,  and  climbing,  as  well  as  high  speed  flight.  However,  inside 
FRA3,  a  factor  Is  introduced  to  make  the  induced  velocity  calculation  more 
flexible  due  to  the  apparent  discrepancy  between  Equations  50  -  51  and  their 
counterparts  which  appear  in  a  report  of  the  Bell  Helicopter  Company  (Refer¬ 
ence  8). 


Whenever  the  table  option  for  Induced  velocity  is  chosen,  an  outside  wake  pro¬ 
gram  has  to  be  run  first  and  the  result  stored  in  a  sequential  file,  the 
appropriate  format  of  which  is  depicted  in  the  User's  Manual  (Volume  II). 


2.3  MODAL  LIFTING  SURFACE  -  CLS2 


A  lifting  surface  structure  is  represented  as  a  linear  combination  of  orthog¬ 
onal  modes.  These  modes  are  defined  as  displacement  and  slope  functions  of 
the  structure  spanwise  coordinate.  The  displacement  at  some  point  on  the  sur¬ 
face  Is  given  by  the  summation  of  the  products  of  the  modal  displacements  at 
that  point  and  the  associated  time-dependent  generalized  coordinates  (modal 
amplitudes): 


l  q,(t) 

i-1  1J  1 


where  Wj  is  the  total  displacement  at  the  point  of  interest,  ^  is  the  ith 
modal  displacement,  and  q^(t)  is  the  ith  generalized  coordinate.  The  slope  is 
given  by 


"}*,!,♦  ij  ‘Mt) 


/  / 

A  control  surface  angle  degree  of  freedom  is  defined  independently  of  the 
modal  degrees  of  freedom  an/ is  coupled  to  them  through  implicit  coupling. 
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/  / 


2.4  LIFTING  SURFACE  AERODYNAMICS  -  FLA2 

The  lifting  surface  aerodynamics  algorithm  is  based  on  the  subsonic  lifting 
surface  theory  of  H.  Multhopp  (References  15  and  16).  This  is  a  theoretical 
method  of  analyzing  finite  wings  that  has  received  wide  acceptance  for  its 
accuracy  and  adaptability  to  a  wide  range  of  fixed-wing  applications  (Refer¬ 
ences  17  and  18). 

The  algorithm  treats  what  is  commonly  referred  to  as  the  inverse  problem: 
i.e.,  given  the  shape  and  attitude  of  a  wing  of  finite  aspect  ratio,  find  the 
corresponding  loading.  A  number  of  features  have  been  included  which  allow 
the  user  much  latitude  in  its  application.  Aerodynamics  forces  can  be  com¬ 
puted  for  full  wings  of  two  semispans  or  for  wings  of  a  single  semispan  (as  in 
the  case  of  a  vertical  tail).  The  wing  may  be  swept  or  unswept,  tapered  or 
untapered,  with  or  without  twist.  The  user  retains  full  control  of  the  choice 
of  the  wing  section  aerodynamic  characteristics  and  is  free  to  vary  them  in 
the  spanwise  direction. 

The  algorithm  Includes  the  option  of  control  surfaces  deflected  symmetrically 
or  antisymmetrically  for  cases  involving  two  semi  spans.  The  effects  of  flaps 
or  ailerons  with  chords  of  approximately  20  percent  of  the  total  sectional 
chord  may  be  computed.  The  user  is  free  to  specify  arbitrarily  the  spanwise 
location  of  the  limits  of  these  control  surfaces. 

The  algorithm  automatically  adjusts  for  the  loss  of  aerodynamic  efficiency  at 
the  spanwise  limits  of  the  control  surface  when  those  spanwise  limits  occur 
within  a  tolerance  of  the  pivotal  stations  at  which  the  downwash  is  computed. 
The  loss  of  aerodynamic  efficiency  is  due  to  the  absence  of  endplating  in 
usual  control  surface  applications.  This  adjustment  takes  the  form  of  a  lin¬ 
ear  reduction  of  the  control  surface  deflection  such  that  if  the  edge  of  the 
control  surface  were  to  coincide  with  the  spanwise  pivotal  point,  the  control 
surface  deflection  would  be  adjusted  to  one-half  its  Inputted  value.  The  span- 
wise  range  over  which  this  adjustment  occurs  varies  from  about  10  percent  of 
the  semispan,  Inboard,  to  about  1.7  percent  at  the  tip. 
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Input  consists  of  the  parameters  describing  the  three-dimensional  wing  plan- 
form  (aspect  ratio,  taper  ratio,  sweep),  two-dimensional  angle  of  attack  and 
sectional  properties  (including  control  surface  characteristics)  at  each  of 
many  spanwise  stations,  plus  the  free-stream  conditions.  The  coordinates  of 
the  wing  are  transformed  by  the  Prandtl -Glauert  compressibility  law  to  allow 
applicability  of  the  solution  up  to  the  critical  Mach  number. 

A  matrix  of  2  chordwise  and  15  spanwise  pivotal  points  is  automatically  com¬ 
puted  over  the  extent  of  the  planform  (see  Figure  6).  It  is  at  these  loca¬ 
tions  that  the  integral  equation  for  the  downwash  is  evaluated.  The  choice  of 
two  chordwise  points  (at  34.55  and  90.45  percent  chord)  at  each  spanwise  sta¬ 
tion  is  necessary  to  fully  incorporate  the  information  of  the  sectional  lift 
and  moment  coefficients  in  the  loading  distribution.  The  choice  of  the  number 
of  spanwise  stations  of  15  is  adequate  for  wings  of  simple  planform  geometry 
and  of  moderate  aspect  ratio.  H.  Multhopp  suggested,  in  Reference  15,  that 
for  two  chordwise  points  the  number  of  spanwise  points,  M,  exceed  the  minimum 
expressed  in  Equation  89  in  terms  of  aspect  ratio,  AR,  and  Mach  number: 


M  >  3  *  AR  *  J\  -  Mach2  (89) 

In  Reference  18,  J.  E.  Lamar  studied  the  effect  of  varying  the  number  of  piv¬ 
otal  points  and  found  that  the  aerodynamic  center  was  the  parameter  most  sen¬ 
sitive  to  the  number  of  these  points.  The  following  formula  was  offered  for 
determining  the  optimal  number  of  pivotal  points  (N  X  M),  where  N  is  the  num¬ 
ber  of  chordwise  points: 


M  -  (4  TO  5)  *  AR  *  J  *  J\  -  Mach2  (90) 

This  equation  (which  does  not  include  sweep)  indicates  that,  for  aspect  ratios 
from  4  to  9  and  Mach  numbers  up  to  0.6,  the  2  X  15  combination  should  be  an 
acceptable  choice  for  most,  but  highly  specialized  planforms,  i.e.,  very  low 
and  very  high  aspect  ratios  with  highly  swept  or  otherwise  specialized  leading 
edge  shapes. 
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STATION  NO. 

NONDIMENSIONAL 

STATION 

1 

0. 

2 

0.1951 

3 

0.3827 

4 

0.5556 

5 

0.7071 

6 

0.8315 

7 

0.9239 

8 

0.9808 

Figure  6.  Lifting  Surface  Pivotal  Points. 


Selection  of  somewhat  less  than  the  optimal  number  of  spanwise  stations  was 
not  found  to  significantly  detract  from  the  results. 

The  lift  and  moment  coefficients  are  automatically  determined  from  the  angle 
of  attack  at  each  of  the  spanwise  stations  computed  by  the  dynamics  package, 
and  from  the  user-provided  aerodynamic  characteristics  (lift  and  moment  linear 
slopes  versus  angle  of  attack,  control  surface  effectiveness,  and  control  sur¬ 
face  deflection). 


The  solution  consists  of  representing  the  mean  cambered  surface  of  the  wing  by 
a  sheet  of  doublets  using  fluid  potential  theory.  The  strengths  of  the  doub¬ 
lets  are  subsequently  computed  at  each  of  the  pivotal  points  of  the  planform 
(as  indicated  in  Figure  6).  Euler's  equation  is  used  to  solve  the  downwash, 
w,  resulting  from  the  loading  distribution,  £(x0,yQ),  where  xQ  and  yQ  are 
points  in  the  planform  as  in  Equation  91: 


w(x,y)  =  gj 


f  ^xp’ 

(y  -  y0>2 


1  + 


(X  -  xn) 


7(x  -  xn)2  +  (1  -  Mach2)(y  -  yn)2 


dx0dy0 

(91) 


With  the  solution  of  the  downwash  field  of  the  wing,  the  induced  angle  of 
attack,  and  then  the  induced  drag,  can  be  computed  by  equations  (92)  and  (93), 
respectively: 


^(x.y) 


w(x,.y) 

U 


(92) 


rb/2 

ci  *  c  •  ai  dy 

>  b/2 


(93) 
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Thus,  the  solution  involves  the  evaluation  of  the  effect  of  each  sectional 
station  on  the  aerodynamic  performance  of  every  other  section  of  the  complete 
planform. 


Special  features  incorporated  in  Reference  15  are  employed  such  as  the  adjust¬ 
ment  in  the  discontinuity  of  the  sweep  at  the  center  line  of  swept  wings,  the 
correction  of  the  logarithmic  singularity  in  the  influence  functions,  the 
choice  of  the  interpolation  polynomials,  etc.  Chordwise  integration  of  the 
two  influence  functions  is  performed  by  a  subroutine  that  uses  Weddle's  rule 
using  24  angular  intervals,  which  was  found  to  be  adequate.  The  solution  in 
the  spanwise  direction  was  carried  out  by  the  iterative  solution  of  successive 
differences  in  the  circulation  and  moment  per  unit  span  as  described  in  Refer¬ 
ence  18.  The  number  of  iterations  was  limited  to  12,  and  the  truncation  terms 
were  then  added  into  the  final  summation  of  circulation  and  moment  per  unit 
span.  From  these  sums,  the  integrated  lift,  drag,  pitching  and  rolling  mom¬ 
ents,  and  center  of  pressure  for  the  given  planform  were  computed.  Inasmuch  as 
the  Multhopp  method  is  well  documented,  more  detailed  discussion  of  the  solu¬ 
tion  of  Equation  91  has  not  been  included. 


In  addition  to  the  limits  of  applicability  noted  above,  the  user  should  not 
expect  to  apply  this  algorithm  to  conditions  that  produce  flow  separation, 
shock  waves  or  otherwise  violate  the  limits  of  linearity  engendered  in  the  use 
of  small  perturbation  analysis. 


2.5  3-D  MODAL  STRUCTURE  -  CFM3 


A  three-dimensional  structure  is  represented  as  a  linear  combination  of  orth¬ 
ogonal  modes.  Modal  translations  and  rotations  are  specified  at  points  de¬ 
fined  in  a  Cartesian  coordinate  system.  The  displacement  at  a  point  in  the 
structure  is  given  by  the  summation  of  the  products  of  the  modal  displacements 
at  that  point  and  the  associated  time-dependent  generalized  coordinates  (modal 
amplitudes) : 
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where  w^  is  the  total  displacement  at  the  point  of  interest,  is  the  ith 
modal  displacement,  and  q.(t)  is  the  ith  generalized  coordinate. 

2.6  LINEAR  CONSTRAINTS  -  CLC2 

Given  a  set  of  arbitrary  relationships  between  degrees  of  freedom, 

[D]{x}-{0),  m<n 
mxn 


where  m  is  the  number  of  constraint  equations  and  n  is  the  number  of  degrees 
of  freedom,  m  implicit  degrees  of  freedom  can  be  solved  for  (eliminated),  m 
degrees  of  freedom  are  selected  from  x,  and  D  is  partitioned  for  retained  and 
eliminated  degrees  of  freedom: 
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Solving  for  Xg: 
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2.7  GENERAL  FORCE  -  CGF2 

All  three  types  of  forcing  function  (see  paragraph  3.1.13,  Volume  II)  are 
applied  periodically.  Start  and  end  times  define  the  period  of  each  force 
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application,  and  the  cycle  is  repeated  automatically  in  the  following  manner. 
Tj  and  T2  are  the  specified  start  and  end  times,  and  T  is  the  current  time. 
The  time  argument  is  defined  as 


T  -  remainder 


-  (T  -  Tj) 


<T2  -  Tl> 


where  the  bracketed  quantity  is  the  largest  integer  whose  magnitude  does  not 
exceed  the  magnitude  of  (T  -  Tj)/(T2  -  Tj). 

2.8  DAMAGED  (NONIDENTICAL)  ROTOR  BLADES  -  CRD3 

The  CRD3  equations  of  motion  are  identical  to  those  of  CRE3.  Changes  in  the 
properties  of  a  blade  or  blades  of  an  existing  rotor  are  specified,  and  the 
coefficients  of  the  component  equations  of  motion  computed  from  the  changes 
are  explicitly  coupled  to  the  original  system. 

Air  loads  are  computed  separately  for  the  damaged  blade(s). 

2.9  LANDING  GEAR  -  CLG2 

A  schematic  of  the  landing  gear  is  shown  in  Figure  7.  x  ,  y  ,  z  ,  a  ,  a  are 

j  j  j  A  jf 

the  degrees  of  freedom  of  the  strut  attachment  point.  z$  is  along  the  axis  of 
the  strut.  The  strut  oleo  is  modeled  as  an  elongation  degree  of  freedom,  AL, 
elastically  coupled  to  the  attachment  degrees  of  freedom,  x^,  y^,  zt  are  tire 
degrees  of  freedom  rigidly  coupled  to  AL.  The  value  of  zt  determines  whether 
or  not  the  tire  is  in  contact  with  the  ground.  When  contact  is  made,  the  tire 
degrees  of  freedom  are  elastically  coupled  to  the  ground  (vertical).  Further, 
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if  friction  is  considered,  the  tire  degrees  of. freedom  are  elastically  coupled 
to  longitudinal  and  lateral  scrubbing  degrees  of  freedom,  x^,  y^.  The  mass  of 
the  strut  is  assigned  to  the  attachment  degrees  of  freedom  and  the  tire  mass 
is  assigned  to  AL  and  the  scrubbing  degrees  of  freedom. 

The  tire  and  scrubbing  degrees  of  freedom  are  defined  with  respect  to  some 
fuselage  coordinate  system.  The  strut  degrees  of  freedom  need  not  be  defined 
in  the  fuselage  coordinate  system  and  must  be  transformed.  (X$,  Y$,  Z$)  are 
the  local  coordinates  of  the  attachment  point.  Writing  the  strut  vectors  in 
terms  of  fuselage  coordinates, 

*S  -  V  +  v  + 

■  bxT  +  byJ  +  b2* 

-*  -t  -t  ■* 

zs  *  V  +  cyJ  +  czk 

where  7,  j,  1c  are  the  unit  vectors  of  the  fuselage  coordinate  system.  Letting 
(xt*  yt,  zt)  represent  fuselage  coordinates,  the  transformation  can  be  written 


and,  accounting  for  translation  and  rotation  of  the  strut, 
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where  L  Is  the  undeformed  length  of  the  strut  (small  displacements  assumed). 
Note  that  xt,  yt,  zt  can  be  treated  as  implicit  degrees  of  freedom. 

\ 

The  transformation  is  then  applied  to  the  component  coefficient  matrices. 

During  a  time  history  solution,  the  transformation  is  applied  to  the  component 
state  vector  to  determine  whether  the  tire  is  in  contact  with  the  ground  and 
to  determine  the  tire  damping  and  stiffness  coefficients  as  a  function  of  ve¬ 
locity  and  displacement.  Strut  damping  and  stiffness  is  computed  directly  from 
the  state  vector. 

When  Coulomb  friction  is  greater  than  or  equal  to  the  stiffness  and  damping 
forces  of  the  tire,  the  tire  is  considered  to  be  elastically  coupled  to  the 
ground  (no  scrubbing);  when  Coulomb  friction  is  less  than  the  combined  stiff¬ 
ness  and  damping  forces,  the  tire  scrubbing  force  is  equal  to  the  friction: 

Lffifliludlnal  i  brates  on 


If:  ux  kz  lztl  *  l°x  (xf  '  xt>  +  kx  (xf  ‘  xt>l 

then:  »  n»t  xf  +  cx  (xf  -  xt)  +  kx  (xf  -  xt) 

else:  fx  *  +  ux  kz  lzt I  ’  xf  <  0 

fx  -  -  ux  kz  lztl  ,  xf  >  0  (u  -  coefficient  of  friction) 

lateral 

If:  uy  kz  |zt|  ^  |cy  (yf  -  yt)  +  ky  (yf  -  yt)| 

then:  f  ■  mt  yf  +  cy  (yf  -  yt)  +  ky  (yf  -  yt) 

else:  f  -  -•-  uy  kz  |zt|  ,  y^  <  0 

f  -  -  u  k  |zf|  ,  yf  ^  0  (u  ®  coefficient  of  friction) 

y  y  z  t  t  y 
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2.10  GENERAL  TINE  HISTORY  SOLUTION  -  STH4 


STH4  performs  a  Runge-Kutta  integration  of  the  system  equations  of  motion. 
Input  parameters  are: 

Start  time 
Initial  increment 
End  time 

Error  check  value 

Initial  velocity,  displacement  of  system  DOF 

Evaluation  is  done  by  means  of  fourth-order  Runge-Kutta  formulae  with  the  mod¬ 
ification  of  Gill  (Reference  19).  Accuracy  is  tested  by  comparing  the  results 
of  the  procedure  for  single  and  double  increments.  The  increment  is  adjusted 
during  the  solution  by  halving  or  doubling.  If  more  than  10  bisections  of  the 
increment  are  necessary  to  achieve  satisfactory  accuracy,  the  solution  is  ter¬ 
minated. 

A  component  state  vector  can  be  obtained  by 
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where  Yj  and  Y  are  the  component  and  system  state  vectors  at  a  given  time  and 
Tj  is  the  component  transformation  matrix. 

2.11  COMPONENT  INTERFACE  AND  INTERNAL  LOADS  -  SI  13  N 

The  complete  equations  of  motion  for  a  component  are 
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where  Mj,  Cj,  and  Kj  are  the  coefficient  matrices  of  the  Ith  component,  Fj  is 
the  forcing  function  vector,  and  FIR  is  the  vector  of  the  reaction  loads  at 
the  interfaces  to  other  components.  All  of  the  FJR  cancel  out  when  a  system 
is  assembled  and  thus  are  not  included  in  the  component  formulation. 

SI 13  computes  the  interface  loads  acting  on  a  component  at  a  given  time: 
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X  and  F  are  the  system  displacement  and  force  vectors  and  Tj  is  the  component 
transformation  matrix. 

Given  Mj,  C j ,  Kj,  Xj,  Xj,  and  Xj,  internal  parameters  can  then  be  computed 
(see  paragraph  3.3.9,  Volume  II). 
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2.12  GENERAL  EIGENANALYSIS  -  SEA5 


The  equation 


[  H  ]  (  X  )  +  [  C  ]  {  X  )  +  [  K  ]  {  X  }  -  {  0  }  . 

nxn 


can  be  transformed  into 


[M](Y}  +  [K]{Y}-{0) 


where 


/  ' 


r  n 


(  Y  } 


x 

X 


[  M  ] 


[  0  ]  [Ml 
CM]  [  C  ] 


[  K  ] 


■CM]  [  0  ] 
[  o  ]  [  K  ] 


The  eigenvalue  problem  is  defined  by 


(Reference  20) 


[  K  ]  {  Y  )  -  [  M  ]  {  Y  }  ,  r  «  1,  2 . 2n 


where  K  and  M  are  real  and  may  be  nonsymmetric,  and  are  the  eigenvalues. 


SEA5  computes  the  eigenvalues  and  eigenvectors  Yr  of  the  above  equation.  K  is 
reduced  to  upper  Hessenberg  form  and  M  to  upper  triangular  form.  K  is  further 
transformed  to  quasi-upper  triangular  form  (upper  Hessenberg  with  no  two  con¬ 
secutive  subdiagonal  elements  being  nonzero).  The  rth  eigenvalue  is  the  com¬ 
plex  number  given  by  the  ratio  of  the  rth  diagonal  elements  of  the  transformed 
K  and  M  matrices  (Reference  21). 
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3.0  INTERFACE  TO  IMPROVED  NONUN ffORM  INFLOW  ANALYSES 

A  dummy  module  called  FRWjJ  has  been  created  to  handle  the  interface  between 
DYSCO  and  a  wake  program.  FRW0  provides  the  following  data:  flapping  angle, 
azimuth  angle  of  blades,  blade  radius,  flapping  hinge  location,  number  of  cal¬ 
culation  points,  nondimensional  radial  location  of  calculation  points,  average 
chord,  average  thickness  ratio,  advance  ratio  (shaft  reference),  inflow  ratio 
(shaft  reference),  air  density,  rotor  speed,  rotor  thrust,  vortex  core  size, 
radial  location  of  root  vortex  (nondimensional),  and  radial  location  of  tip 
vortex  (nondimensional). 

Like  other  aerodynamic  force  modules,  i.e.,  FRA0,  FRA2  and  FRA3,  FRW0  may  be 
coupled  with  the  rotor  modules  CRR2  and  CRE3  which,  in  turn,  may  be  coupled 
with  any  other  component  modules  to  form  a  model  within  the  limitations  of 
DYSCO. 

FRW0  consists  of  four  modules:  input  module  FRW0I,  constant  module  FRW0C, 
active  module  FRW0A,  and  module  FRW0B.  A  more  detailed  description  of  each  of 
these  modules  is  given  in  the  following  paragraphs. 

3.1  HQPUIE.  FRW0I 

The  values  of  the  following  parameters  in  this  module  must  be  provided  by  the 
user: 

1.  Component  of  wind  velocity,  with  regard  to  hub  system  (for  a 
trim  case,  the  wind  velocity  should  be  set  to  zero  because  it 
Is  replaced  by  vehicle  velocity  and  wind  velocity  with  regard 
to  ground  system  in  the  Input  module  of  trim  solution  input 
module) 

2.  Number  of  nondimensional  stations  (calculation  points) 

3.  Nondimensional  stations  (r/R) 

4..  Rotor  thrust 

i 


83 


5.  Average  thickness  ratio 

6.  Vortex  core  size 

7.  Nondimensional  station  of  root  vortex 

8.  Nondimensional  station  of  tip  vortex. 

3.2  MODULF  FRW0C  ,  , 

/ 

In  this  module,  all  input  data  are  retrieved  and  stored.  Information  con¬ 
cerning  wind  velocity  and  rotor  thrust,  air  density  ratio,  average  thickness 
ratio,  vortex  core  size,  station  of  root  vortex  and  tip  vortex  are  stored  in 
common  /FRW0/,  and  data  concerning  number  of  stations  and  nondimensional  sta¬ 
tions  are  stored  in  common  /ROT/.  Both  /FRW0/  and  /ROT/  also  appear  in  module 
FRW0A. 

If  coupled  with  elastic  rotor  CRE3,  this  module  also  constructs  mode  shape  and 
blade  pitch  angle  based  on  the  aerodynamic  stations  using  the  interpolation 
method. 

3.3  MODULE  FRW0A 

Using  the  data  in  /FRW0/  and  /ROT/,  calculations  of  advance  ratio  and  inflow 
ratio  about  the  shaft  reference  axis  are  performed  in  this  module.  These  cal¬ 
culations,  together  with  information  about  number  of  blades,  azimuth  angle  of 
blades,  blade  radius,  flapping  hinge  location,  and  flapping  angle  of  each 
blade,  which  according  to  the  DYSCO  convention  should  be  provided  by  the  rotor 
component,  provide  all  the  data  required  by  this  contract. 

3.4  MODULE  FRW0B 

When  the  'CASE'  command  is  activated,  FRW0B  is  called  to  store  all  the  infor¬ 
mation  concerning  /FRW0/  into  a  user  data  file  for  batch  purposes. 
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3.5  SUMMARY 


Details  of  the  computed  parameters  are  listed  below: 


Geometric  Parameters 

NB(NROT)  -  Number  of  blade  of  NROTth  rotor,  provided  by 

rotor  component  stored  in  /ROT/. 

SLOZ(K,L)  -  Flapping  angle  of  Kth  Station  of  Lth  blade, 

provided  by  rotor  component  through  argument 
list. 


PSIN 

R 


EFC(IREF) 

NXA 


X(I) 


CHORDC(I) 

TRTO(FREF) 


-  Azimuth  angle  of  Lth  blade,  provided  by  rotor 
component  through  argument  list. 

-  Blade  radius,  provided  by  rotor  component, 
stored  in  /ROT/  as  XC(NX,  NROT)  where  NX  is 
number  of  dynamic  stations. 

-  Flapping  hinge  location,  provided  by  rotor 
component,  stored  in  /ROT/. 

-  Number  of  aerodynamic  stations  (calculation 
points),  provided  by  FRW0I  stored  in  /ROT/  as 
NAERO(FREF)  where  FREF  is  the  FREFth  use  of 
FRW0. 

-  Radial  location  of  calculation  points,  pro¬ 
vided  by  FRW0I,  stored  in  /ROT/  as 
XAER0(I,NR0T) . 

-  Average  blade  chord,  provided  by  FRW0I, 
stored  in  /ROT/. 

-  Average  thickness  ratio,  provided  by  FRW0I, 
stored  in  /ROT/. 
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Flight  Condition  Parameters 


AMUN 

-  Advance  ratio  (shaft  reference),  calculated 
in  FRWjJA. 

AINFWN 

-  Inflow  ratio  (shaft  reference),  calculated  in 
FRWjJA. 

FACT 

-  Air  density,  provided  by  FRWjJI  and  stored  in 
/FRWjJ/. 

OM(NROT) 

-  Rotor  rpm,  provided  by  rotor  component  and 
stored  in  /ROT/. 

RCT(FREF) 

-  Rotor  thrust,  provided  by  FRWjJI  and  stored  in 
/FRWjJ/  as  RCTC(FREF).  Note  that  for  trim 
case  this  value  should  be  replaced  by  CT  in 
/TRIM/  after  1st  iteration. 

/  V 

Wake  Parameters 

VCS(FREF) 

-  Vortex  core  size,  provided  by  FRWjJ  I  and 
stored  in  /FRWjJ/. 

RVR(FREF) 

-  Radial  location  (nondimensional )  of  root  vor¬ 
tex,  provided  by  FRWjJI  and  stored  in  /FRWjJ/. 

RVT(FREF) 

-  Radial  location  (nondimensional)  of  tip  vor¬ 
tex,  provided  by  FRWjJI  and  stored  in  /FRWjJ/. 
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The  primary  goal  of  this  validation  effort  is  to  compare  DYSCO  computed  trim 
control  parameters  and  horsepower  consumption  with  those  measured  in  the  Oper¬ 
ational  Loads  Survey  test  flight  of  the  Bell  AH-1G  in  order  to  determine  the 
program's  capabilities  and  accuracy  (Reference  22).  In  addition,  the  results 
obtained  from  the  elastic  blade  analysis,  using  both  the  advanced  aerodynamic 
analysis  and  the  simplified  aerodynamic  analyses  currently  in  DYSCO,  are  com¬ 
pared  with  each  other  to  observe  differences  between  them.  The  model  vehicle 
is  the  Bell  AH-1G,  a  teetering  single-rotor  helicopter  with  anti-torque  tail 
rotor. 

The  general  input  data  for  the  AH-1G  associated  with  C81  (Reference  23)  are 
presented  in  the  appendix. 

4.1  DYSCO  TRIM  ALGORITHM 

In  DYSCO,  the  trim  condition  is  defined  for  a  fuselage  flying  at  a  steady 
state  condition  with  the  blades  and  fuselage  in  periodic  motion.  The  trim 
algorithm  by  D.  A.  Peters  and  A.  P.  Izadpanah,  called  "Periodic  Shooting  With 
Newton-Raphson  Iteration,"  Reference  24,  has  been  adapted  to  include  a  com¬ 
plete  helicopter,  nonlinear  aerodynamics,  and  such  flight  conditions  as  steady 
turn,  steady  climb,  and  autorotation.  This  method  requires  only  numerical 
integration  of  nonlinear  equations  and  the  inversion  of  an  approximate 
partial -derivative  matrix  formed  from  successive  numerical  integrations 
through  one  period  for  every  perturbation  in  controls  and  initial  conditions. 
This  algorithm  has  the  following  characteristics: 

1.  On  a  computational  basis,  it  is  superior  to  any  other  transi¬ 
tion-matrix  method  in  that  it  requires  less  computer  storage 
and  computation. 

2.  It  is  superior  to  the  harmonic  balance  method,  except  for  a 
linear  system  with  few  harmonics  in  the  coefficients. 

3.  It  is  superior  to  a  direct  numerical  integration  for  a  case 
with  small  damping. 
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As  far  as  control  strategy  is  concerned,  four  independent  direct  controls  and 
two  indirect  controls  are  used.  The  choice  of  direct  controls  depends  on  the 
rotor  type  and  flight  condition.  The  two  indirect  controls  are  two  of  the 
three  Euler  angles  of  a  ground  fixed  coordinate  system  with  respect  to  the 
body  fixed  system.  This  algorithm  treats  all  the  controls  and  initial  condi¬ 
tions  of  blade  degrees  of  freedom  as  unknowns  and  adjusts  them  simultaneously 
to  achieve  force  and  moment  balance  in  the  concerned  directions,  as  well  as 
steady  vibration  of  all  the  degrees  of  freedom  of  the  system. 

The  helicopter  trim  solution  module  (STR3)  involves  an  analysis  which  is  some¬ 
what  more  specialized  than  the  general  dynamic  or  aerodynamic  analyses  often 
carried  out  with  DYSCO.  That  is,  the  trim  algorithm  includes  large  angular 
motions  with  respect  to  a  fixed  coordinate  system,  variations  in  the  relation¬ 
ship  of  the  body  axis  to  the  wind  vector,  and  automatic  treatment  of  the  grav¬ 
ity  vectors.  For  these  reasons,  the  present  applications  require  adherence  to 
certain  considerations  not  necessary  in  other  DYSCO  analyses. 

At  present,  it  is  recommended  that  users  apply  STR3  to  single  dynamic  rotor, 
single  fuselage  combinations  because  a  set  of  independent  trim  variables  must 
be  introduced  for  more  complex  systems. 

42  GENERAL  CONSIDERATION  OF  THE  MODFL  FOR  THE  TRIM  SOLUTION 

4-2-l  CRR2,  CRE3  (Rotor  Component! .  Use  a  single  rotor  only.  The  number  of 
blades  is  arbitrary.  A  teetering  rotor  may  be  modeled  by  the  additional  use 
of  a  CLC1  module  (BETA1200  -  -1  *  BETA1100  for  CRR2,  or  OP  1110  =»  -1  * 
OP  i210  for  CRE3) .  A  "single  blade  analysis'  may  be  performed  for  efficiency 
of  computations.  The  blade  degrees  of  freedom  are  arbitrary,  except  for  the 
flap  or  out-of-plane  degree  of  freedom.  If  pitch  or  torsion  is  a  DOF,  CCE0  or 
CCE1  must  also  be  used  to  provide  controls  for  the  rotor. 

All  other  parameters  are  optional;  the  hub  degrees  of  freedom  must  be  consis¬ 
tent  with  the  overall  model. 
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4.2.2  FRAI*  _FRA2^_FRA3  (Rotor  Aerodynamic  Hodulel .  The  wind  velocity  should 
be  entered  as  zero  since  It  Is  an  Input  to  STR3.  Care  must  be  exercised  in 
selecting  the  nondlmensional  stations  for  the  aerodynamics  which  represent  tip 
and  root  loss  effects.  The  safest  procedure,  at  present,  is  to  select  values 
that  correspond  to  the  aerodynamic  stations. 

4.2.3  CFH2  f Fusel aae  Component).  The  first  CFM2  in  the  model  is  the  "refer¬ 
ence"  fuselage.  (At  present,  it  is  recommended  that  only  one  be  used  in  a 
trim  analysis.  Contact  Kaman  for  guidance  If  it  is  desired  to  use  additional 
CFM2  modules.)  The  rotor  should  be  coupled  to  this  component  using  the  stan¬ 
dard  input. 

The  reference  axes  are  the  principal  axes  of  the  fuselage.  Note  that  the  CG 
is  that  of  the  component,  not  of  the  entire  vehicle.  The  mass  and  gravita¬ 
tional  effects  of  the  rotor  and  fuselage  are  automatically  managed  by  STR3. 

The  rigid  body  degrees  of  freedom  selected  determine  which  forces  and  moments 
are  to  be  balanced  during  trim.  Only  three  options  are  presently  allowed: 

1.  XCG,  ZCG,  ROLL,  PTCH 

2.  XCG,  YCG,  ZCG,  ROLL,  PTCH 

3.  XCG,  YCG,  ZCG,  ROLL,  PTCH,  YAW. 

All  options  allow  level  flight  or  steady  climb  or  descent.  Options  2  and  3 
may  also  be  used  for  steady  turn.  Option  3  may  be  used  for  autorotation. 

4.2.4  FFC2  (Fuselage  Aerodynamic  Module).  This  aerodynamic  force  module 
should  be  used  with  the  reference  fuselage.  Wind  velocity  should  be  input  as 
zero,  since  this  is  Input  in  STR3.  When  tall  rotor  thrust  is  a  trim  variable, 
it  is  iterated  upon  and  the  input  value  in  this  module  may  be  an  estimated 
value.  Note  that  when  a  teetering  rotor  is  modeled,  the  horizontal  tail  inci¬ 
dence  is  coupled  to  the  cyclic  rotor  control  (see  STR3  input).  Also  note  that 
all  the  moment  coefficients  in  this  module  are  based  on  the  characteristic 
length  of  1  inch. 
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4*2.5  CLC1  (Linear  Constnlnt).  This  module  provides  the  necessary  con- 
straint  of  the  rigfB  body  flapping  modes  of  the  two  blades  to  form  the  tee¬ 
tering  mode.  In  some  cases,  this  module  may  also  be  used  to  scale  the  elastic 
mode  which  Is  equivalent  to  changing  the  tip  deflection  of  the  mode  shape  to 
another  constant  instead  of  one.  This  may  avoid  numerical  instability  in  some 
situations. 

4*2.6  CCEf.«_.  CC£1.  (Rotor  Control  Svsteml.  The  use  of  a  control  system  is  man¬ 
datory  whenever  the  rotor  pitch  or  elastic  twist  degrees  of  freedom  are  se¬ 
lected. 


4.2.7  Other  Components.  Any  other  components  having  a  negligible  mass  may 
also  be  Included  In  the  model. 

4.2.8  STR3  Input.  The  following  comments  should  assist  the  user  In  supplying 
appropriate  Input  values: 

If  teetering  rotor  -  Input  A  and  8  such  that  the  horizontal  tail 
Incidence  angle  -  A  +  B  *  cyclic  sine  control  +  (C  *  cyclic  sine 
control )**2  (radians). 

Initial  Integration  Increment  (seel  -  Depends  on  system  frequencies. 

If  only  rigid  body  modes,  an  appropriate  value  Is  1/20  of  the  time 
for  one  rotor  revolution. 

Integration  Period  (seel  -  Normally  the  time  for  one  rotor  revolu¬ 
tion. 

Mfl.t  of.  Iterations  -  5  -  10  is  typical  . 

Error  Allowed  for  Each  DOF  -  i.e.,  allowable  tolerance  (absolute). 

Increment  on  Each  DOF  -  For  calculating  partial  derivative  matrix. 

.01  radian  or  .01  Inch  are  typical. 

Wind  Velocity  -  In  terms  of  ground  fixed  system  (In./sec).  Absolute 
velocity  of  air  with  respect  to  the  ground  (not  with  respect  to  the 
fuselage). 
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Fuselage  C6  Translational  Velocity  tin, /sec I  -  In  terms  of  ground 
fixed  system.  Include  tangential  velocity  for  a  turn.  Note  X  velo¬ 
city  Is  negative  In  forward  flight  (see  CFM2  coordinate  system  in 
User's  Manual). 

Fuselage  Rotational  Velocity  (rad/sec)  -  Absolute  rotational  velo¬ 
city  of  fuselage  In  terms  of  ground  fixed  system. 

Initial  Estimates  -  May  be  zero  for  initial  conditions  and  controls. 

4.3  TRIM  SIMULATION  OF  AH-1G 

As  mentioned  in  the  previous  section,  the  DYSCO  trim  algorithm  involves  an 
analysis  which  is  somewhat  more  specialized  than  the  general  dynamic  or  aero¬ 
dynamic  analyses  often  carried  out  with  DYSCO.  So,  when  a  trim  model  is 
formed,  one  has  to  follow  the  instructions  of  the  previous  section  or  the 
User's  Manual  carefully  in  order  to  avoid  any  unnecessary  mistakes. 

The  simulation  of  the  AH-1G  by  DYSCO  Includes  several  component  and  force 
modules:  elastic  rotor  CRE3,  rotor  aerodynamic  module  FRA3,  rotor  control 
system  CCE0,  linear  constraint  module  CLC1,  fuselage  component  CFM2,  and  fuse¬ 
lage  aerodynamic  module  FFC2. 

The  following  Is  a  detailed  description  of  the  simulation  and  its  limitations. 

Because  the  derivation  of  the  equations  of  motion  in  the  aeroelastic  rotor 
component  CRE3  Is  based  on  the  energy  method  (or  its  equivalent,  the  general¬ 
ized  Galerkin  method),  the  Input  modes  need  not  be  orthogonal  as  long  as  they 
satisfy  the  geometric  boundary  conditions.  In  addition  to  this,  these  modes 
may  be  obtained  for  each  physical  degree  of  freedom  separately.  In  this  simu¬ 
lation,  the  In-plane  and  out-of-plane  elastic  modes  were  obtained  from  the 
eigensolutlon  of  a  B540  blade  rotating  at  324  rpm  with  fixed-free  boundary 
conditions.  The  torsional  modes  consist  of  the  rigid  body  mode  and  a  linearly 
varied  elastic  mode.  All  of  these  modes  are  superimposed  on  the  out-of-plane 
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rigid  body  lode  which  consists  of  the  inodes  necessary  for  a  teetering  rotor. 
Note' that  it  is  thfc  space  spanned  by  these  modes,  and  not  the  Individual  mode, 
that  dominates  the  accuracy  of  the  result. 

Although  all  data  concerning  B540  blades  are  available  in  Reference  23,  the 
rectangular  input  for  blade  distributed  properties  in  Reference  23  Is  diffi¬ 
cult  to  convert  to  the  data  needed  for  OYSCO,  which  accepts  only  the  trape¬ 
zoidal  distributed  properties. 

One  way  to  avoid  this  problem  Is  to  Introduce  double  stations  for  each  radial 
location.  However,  this  may  greatly  affect  computational  efficiency  and  is 
not  adopted  here.  Zero  value  for  the  precone  angle  is  used  in  this  simula¬ 
tion  because  the  current  version  of  DYSCO  cannot  handle  the  transformation 
when  two  degrees  of  freedom  are  differentiated  only  by  a  constant.  However, 
it  is  expected  that  the  teetering  constraint  of  the  two  blades  may  totally 
cancel  the  effect  of  the  precone  angle  to  the  linear  terms  and  have  negligible 
effect  on  the  trim  result. 

The  simplest  control  system  CCE0  which  considers  only  the  control  rod  stiff¬ 
ness  is  used  to  simulate  the  AH-1G  control  system  because  it  Is  a  recognized 
fact  that  an  equivalent  spring  may  be  used  to  reasonably  represent  a  helicop¬ 
ter  control  system. 

In  CFM2,  four  rigid  body  modes  (XCG,  ZCG,  ROLL,  and  PTCH)  are  chosen  and  no 
elastic  mode  is  considered. 

DYSCO  can  perform  a  trim  solution  for  the  three  force  balances  and  the  three 
moment  balances  of  the  fuselage.  However,  because  only  the  two-dimensional 
wind  effect  Is  considered  for  the  aerodynamic  surfaces  In  the  fuselage  aerody¬ 
namic  force  module  FFC2  and  because  a  tail  rotor  is  not  allowed  by  the  trim 
algorithm,  the  Y-dlrectlon  force  balance  and  yaw  moment  balance  of  the  fuse¬ 
lage  are  not  considered  In  this  simulation.  This  implies  that  only  four 
degrees  of  freedom  -  XCG,  ZCG,  PTCH,  and  ROLL  -  are  assigned  to  the  fuselage. 

\  \ 
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In  CFM2,  the  rotor  shaft  angles  (forward  and  lateral)  from  vertical  and  the 
vertical  height  and  horizontal  distance  of  the  hub  from  the  CG  provide  all  of 
the  information  needed  to  automatically  couple  with  rotors  through  rigid  mass¬ 
less  shafts. 

Note  that  In  DYSCO  the  fuselage  CG  does  not  include  the  contribution  from  the 
rotor  and  Is  different  from  the  vehicle  CG,  which  is  used  for  helicopter  oper¬ 
ation. 

Note  also  that  the  mass  moments  of  inertia  of  the  fuselage,  which  are  required 
as  input,  are  not  used  in  trim. 

A  teetering  rotor  can  be  formed  by  using  module  CLC1  to  constrain  the  rigid 
flapping  modes  of  the  two  blades.  This  module  can  also  be  used  to  scale  the 
coefficient  matrices  and  make  the  numerical  procedure  more  stable. 

Fuselage  aerodynamic  force  module  FFC2  only  considers  the  two-dimensional 
aerodynamic  effect  by  expressing  each  aerodynamic  coefficient  as  a  quadratic 
function  of  angle  of  attack  In  the  plane  of  the  airfoil  cross  section. 

The  aerodynamic  force  on  the  wing  can  be  calculated  using  the  equations  in 
Reference  24,  but  this  method  is  not  yet  available  in  FFC2  for  the  current 
version  of  DYSCO.  A  regression  procedure  is  employed  to  determine  the  wing 
aerodynamic  coefficients  based  on  the  trimmed  fuselage  aerosurface  forces  and 
moments  from  C81  at  various  speeds  of  Interest. 

4.4  RESULTS  AND  DISCUSSION 

The  basic  DYSCO  model  created  in  this  validation  effort  is  AH-1G-35A,  which  is 
based  on  the  AH-1G  helicopter  with  a  gross  weight  of  8300  lb,  aft  CG,  and 
counter  number  35A  used  in  the  operational  loads  survey  test  flight  (Refer¬ 
ences  22  and  23).  The  general  characteristics  of  the  AH-1G  helicopter  are 
presented  in  the  appendix.  The  DYSCO  trim  model  AH-1G-35A  is  demonstrated  in 
Figure  8. 


*************  MODEL  AH-1G-35A  ************** 
AH-1G  TRIM  SAMPLE  CASE 


INDEX 

COMP 

NO. 

DATA  SET 

FORCE 

DATA  SET 

1 

CRE3 

1 

B2Z1T2 

FRA3 

FCT1.65 

REQUIRED  DS/DM-AFD161/AIRI 

2 

CCEO 

1 

3000 

NONE 

3 

CLC1 

COUPLE 

NONE 

4 

CFM2 

1 

8300-4 

FFC2 

AH1G16.5 

*********************************** 


Figure  8.  AH- 1G-35A  Trim  Model . 


Detailed  Information  for  component  and  force  modules  associated  with  the 
AH-1G-35A  model  can  be  seen  in  the  User's  Manual.  However,  for  the  sake  of 
clarity,  some  general  descriptions  of  these  components  and  forces  are  listed 
below: 

1.  CRE3,  B1Z1T2  -  2  blades;  In-plane  DOFs,  one  mode  (IP  1110, 

IP  1210);  out-of-plane  DOFs,  two  modes  (OP  1110,  OP  1120, 

OP  1210,  OP  1220);  torsion  DOFs,  two  modes  (TOR  1110, 

TOR  1120,  TOR  1210,  TOR  1220);  4  hub  DOFs  (XHUB1000,  ZHUB1000, 
ALFX1000,  ALFY1000);  14  degrees  of  freedom  total. 

2.  FRA3,  FCT1.65  -  Induced  velocity  by  equation  with  a  1.65  fac¬ 
tor,  aerodynamic  coefficients  by  table  look-up.  Unsteady  Is 
off. 

3.  CCE0,  3000  -  Control  system  with  3000  lb/in.  control  rod  stiff¬ 
ness,  coupled  with  blade  pitch  or  torsion  degrees  of  freedom. 
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4.  CLC1,  COUPLE  -  Teetering  constraint  for  out-of-plane  mode 
(OP  1110»264*TEET,  OP  1210— 264*TEET) ;  scaling  for  first 
in-plane  mode  and  second  out-of-plane  mode 

(IP  1110»264*IPIP1110,  IP  1210«264*IPIP1210, 

OP  1 1 20»264*OPOP1 1 20 ,  OP  1220-264*OPOP1220) . 

5.  CFM2,  8300-4  -  4  rigid  body  modes  (XCG  1000,  ZCG  1000, 

ROLLIOOO,  PTCH1000);  automatically  couples  to  rotor; 

mass  «  7290  lb,  station  of  hub  ■  -  .68  in.,  vertical  height  of 
hub  -  96.485  In. 

6.  FFC2,  AH1G16.5  -  Flat  plate  drag  area  is  16.5  square  ft.  with 
fuselage,  wing,  horizontal  tail,  and  vertical  tail  aerody¬ 
namics. 

Another  DYSCO  model  AH-1G-36A  corresponding  to  a  9000-lb  gross  weight,  mid-CG 
location,  and  counter  number  36A  in  the  operational  loads  survey  test  flight 
is  shown  in  Figure  9. 

*************  MODEL  AH-1G-36A  ************** 

AH-1G  TRIM  SAMPLE  CASE 


INDEX 

COMP 

NO. 

DATA  SET 

FORCE 

DATA  SET 

1 

CRE3 

1 

B2Z1T2 

FRA3 

PCT1.65 

REQUIRED  DS/DM-AFD161/AIRF0IL 

2 

CCEO 

1 

3000 

NONE 

3 

CLC1 

COUPLE 

NONE 

4 

CFM2 

1 

9000-A 

FFC2 

AH1G16.5 

*********************************** 


Figure  9  .  AH-1G-36A  Trim  Model. 
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Model  AH-1G-36A  Is  constructed  by  simply  replacing  8300-4/CFM2  in  AH-1G-35A 
with  9000-4/CFM2,*%hich  is  the  same  as  8300-4/CFM2  except  that  mass  -  7990  lb, 
station  of  hub  -  4.278  in.,  and  vertical  height  of  hub  -  95.427  in.  Further¬ 
more,  in  order  to  compare  the  differences  between  rigid  rotor  CRR2  and  elastic 
rotor  CRE3  (with  only  rigid  blade  modes)  and  among  aerodynamic  modules  FRA0, 
FRA2  and  FRA3,  four  additional  models  were  created  and  are  illustrated  in  Fig¬ 
ures  10  through  13. 

The  component  and  force  data  sets  for  the  four  additional  models  are  described 
below: 

a.  CRE3,  B1Z1T1  -  2  blades;  in-plane  DOFs*  one  rigid  mode 
(IP  1110,  IP  1210);  out-of-plane  DOFs,  one  rigid  mode 
(OP  1110,  OP  1210);  torsion  DOFs,  one  rigid  mode  (TOR  1110, 

TOR  1210);  4  hub  DOFs  (XHUB1000,  ZHUB1000,  ALFX1000,  ALFY1000); 

10  degrees  of  freedom  total. 

b.  CRR2,  BZT  -  2  blades;  lag  DOFs  (ZETA1100,  ZETA1200);  pitch  DOFs 
(THET1100,  THET1200) ;  flapping  DOFs  (BETA1100,  BETA1200);  4  hub 
DOFs  (XHUB1000,  ZHUB1000,  ALFX1000,  ALFY1000),  10  degrees  of 
freedom  total;  other  characteristics  are  the  same  as 
B1Z1T1/CRE3. 

c.  CLC1,  TEET  -  Teetering  constraint  of  rigid  rotor  (TEET  * 

BETA1100  -  -BETA1200). 

d.  CLC1,  TEETE  -  Teetering  constraint  of  elastic  rotor  (TEET  - 

0P1 110/264  0P1 210/264. 

e.  CLC1,  TSCALE  -  Scaling  component  for  elastic  rotor  in-plane 

mode.  (IP  1110  -  264*IPSC1110,  IP  1210  -  264*IPSC1210). 

f.  FRA3,  General  -  Induced  velocity  factor  -  2,  Tiplos  coeffi¬ 

cient  -  1,  other  characteristics  are  the  same  as  FCC1.65/FRA3. 

g.  FRA2,  Tabular  -  Characteristics  are  the  same  as  General/FRA3, 
but  use  simple  momentum  theory  to  calculate  induced  velocity; 
no  yaw  flow  correction  for  aerodynamic  calculation. 
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**************  MODEL  ELAST3  ************** 
ELAST3:  ELASTIC  TEETERING  ROTOR  AND  GENERAL  AERODYNAMICS 


INDEX 

COMP 

NO.  DATA  SET 

FORCE 

DATA  SET 

1 

CRE3 

1  B2Z1T2 

FRA3 

GENERAL 

REQUIRED  DS/DM*AFD161/AIRF0IL 

2 

CCEO 

1  3000 

NONE 

3 

CLC1 

TEETE 

NONE 

4 

CLC1 

TSCALE 

NONE 

5 

CFM2 

1  8300-4 

FFC2 

AH1G16.5 

*  *  *  * 

*  *  *  * 

*  *  *.*  *****  *  * 
*\ 

******* 

******** 

Figure  10.  ELAST3  Trim  Model. 

*  *  *  * 

********** 

MODEL  ELAST2 

************** 

ELAST2 : 

ELASTIC  TEETERING 

ROTOR  AND  TABULAR  AERODYNAMICS 

INDEX 

COMP  NO. 

DATA  SET 

FORCE  DATA  SET 

1 

CRE3  1 

B2Z1T1 

FRA3  TABULAR 

REQUIRED  DS/DM-AFD161/AIRF0IL 

2 

CCEO  1 

3000 

NONE 

3 

CLC1 

TEETE 

NONE 

4 

CLC1 

TSCALE 

NONE 

5 

CFM2  1 

8300-4 

FFC2  AH1G16.5 

*****************  ***************** 

Figure  11.  ELAST2  Trim  Model. 
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**************  moq£L  ELASTO  ************** 

ELASTO:  ELASTIC  TEETERING  ROTOR  AND  LINEAR  AERODYNAMICS 


INDEX 

COMP 

NO. 

DATA  SET 

FORCE 

DATA  SET 

1 

CRE3 

1 

B1Z1T1 

FRAO  • 

LINEAR 

2 

CCEO 

1 

3000 

NONE 

3 

CLC1 

TEETE 

NONE 

4 

CLC1 

TSCALE 

NONE 

X  / 

5 

CFM2 

1 

8300-4 

FFC2 

f 

AH1G16.5 

*  *  *  * 

*  *  *  * 

*  *  *  * 

*****  *  *  * 

****** 

****** 

Figure  12.  ELASTO  Trim  Model. 


**************  MODEL  RIGIT2  ************** 
RIGIT2:  RIGID  TEETERING  ROTOR  AND  TABULAR  AERODYNAMICS 

INDEX  COMP  NO.  DATA  SET  FORCE  DATA  SET 

1  CRR2  1  BZT  FRA2  TABULAR 

REQUIRED  DS/DM*AFD1 6 1/A I RFO I L 

2  CCEO  1  3000  NONE 

3  CLC1  TEET  NONE 

4  CLC1  ORDER  NONE 

5  CFM2  1  8300-4  FFC2  AH1G16.5 

*****************  ***************** 

Figure  13.  RIGIT2  Trim  Model . 


98 


h.  FRA0,  Linear  -  Lift  coefficient  -  .112/deg,  drag  coeffi¬ 
cient  »  .085,  moment  coefficient  «  0,  other  characteristics  are 
the  same  as  Tabular/FRA2. 

For  the  AH-1G-35A  case  (corresponding  to  flight  *15A,  clean  wing,  8300  lb  gross 
weight,  aft  CG  station  200.6),  the  comparisons  cf  DYSCO  computed  and  OLS  mea¬ 
sured  flight  performance  parameters  are  shown  in  Figures  14a  through  14e.  In 
these  figures,  the  triangle,  AH-1G-35A(S),  represents  the  DYSCO  model  AH-1G- 
35A  with  steady  aerodynamics;  the  diamond,  AH-1G-35A(U),  represents  DYSCO 
model  AH-1G-35A  with  unsteady  aerodynamics;  and  the  circle,  AH-1G-35A(B1S) , 
stands  for  the  same  DYSCO  model  considering  only  the  rigid  out-of-plane  mode 
with  steady  aerodynamics.  As  can  be  seen  in  these  figures,  the  general  trend 
and  agreement  between  measured  and  computed  results  are  reasonably  good  if  the 
irregular  test  points  are  neglected.  Also,  it  is  observed  that  the  model 
AH-1G-35A(B1$),  which  considers  only  one  rigid  flapping  mode,  also  gives  good 
prediction  as  far  as  trim  positions  are  concerned. 

The  consideration  of  unsteady  aerodynamics  has  virtually  no  effect  on  the 
horsepower  and  has  little  effect  on  fuselage  pitch  or  collective  and  cyclic 
control  angles,  as  demonstrated  in  Figures  14a  through  14e. 

The  effect  of  different  modes  on  the  trim  result  is  depicted  in  Figures  15a 
through  15e.  In  Figures  15,  the  circle,  AH-1G-35A(S),  represents  the  DYSCO 
model  AH-1G-35A  with  steady  aerodynamics  (shown  in  Figure  8);  the  square, 
AH-1G-35A(B1S),  the  triangle,  AH-1G-35A(B35),  and  the  diamond,  AH-1G- 
35A( B3Z25) ,  represent  the  same  model  with  one  rigid  out-of-plane  mode,  three 
out-of-plane  modes,  and  three  out-of-plane  and  two  in-plane  modes,  respec¬ 
tively.  It  is  seen  in  Figure  15a  that  the  collective  control  position  of 
AH-1G-35A(S)  is  displaced  by  a  constant  value  as  compared  with  the  other  three 
models.  The  reason  for  this  is  that  there  is  no  torsional  degree  of  freedom 
for  the  other  three  models;  therefore,  the  control  is  directly  applied  to  the 
blade  which  is  equivalent  to  being  controlled  by  a  control  system  with  infi¬ 
nite  spring  rate.  It  is  illustrated  in  Figures  15a  through  15e  that  the 
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FLIGHT  TEST  35A 
AH1G-35A  (S) 


Comparison  Between  Flight  35A  and  DYSCO  Simulation 
Collective  Control  vs  Advance  Ratio. 
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ADVANCE  RATIO 


Figure  14c.  Comparison  Between  Flight  35A  and  DYSCO  Simulation, 
Cyclic  Sine  Control  vs  Advance  Ratio. 
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AH1G-35A  (BIS) 


Degree  of  Freedom  Effect  on  Trim  -  Flight  35A  Simulation 
Collective  Control  vs  Advance  Ratio. 
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Figure  15b.  Oegree  of  Freedom  Effect  on  Trim  -  Flight  35A  Simulation 
Cyclic  Cosine  Control  vs  Advance  Ratio. 
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low  Effect  on  Trim  -  Flight  35A  Simulation 
Angle  vs  Advance  Ratio 


In-plane  inodes  have  negligible  effect  on  the  rotor  performance  (compare  tri¬ 
angle  and  diamond).  On  the  other  hand,  the  two  elastic  out-of-plane  bending 
modes  have  moderate  effect  on  the  trim  result  (compare  square  and  triangle). 
Despite  the  differences  in  these  four  models  in  predicting  the  trim  position, 
it  is  interesting  to  see  that  they  require  virtually  the  same  amount  of  horse¬ 
power  as  shown  in  Figure  15e.  The  difference  in  cyclic  cosine  control  angles 
for  various  models  as  shown  in  Figure  15b,  however,  requires  further  study. 

The  comparison  between  measured  and  computed  results  associated  with  flight 
36A  is  shown  in  Figures  16a  through  16e.  The  meanings  of  the  triangle,  dia¬ 
mond,  and  circle  in  these  figures  are  the  same  as  those  in  the  AH-1G-35A  case. 
As  can  be  seen  in  Figures  16a,  16c,  16d,  and  16e,  the  agreement  among  measured 
and  computed  results  is  generally  good.  Again,  the  difference  in  cyclic 
cosine  angles,  as  shown  in  Figure  16b,  seems  to  require  further  Investigation. 

Four  DYSCO  models  have  been  created,  as  shown  in  Figures  10  through  13.  In 
the  model  ELAST2,  the  elastic  rotor  component  B1Z1T1  is  assigned  three  rigid 
body  modes,  one  for  each  of  the  three  blade  physical  degrees  of  freedom. 
Therefore,  ELAST2  can  be  used  for  comparison  with  the  results  of  model  RIGIT2 
to  achieve  a  limited  validation. 

It  can  also  be  seen  in  these  figures  that  the  three  aeroelastic  rotor  models 
ELASTO,  ELAST2,  and  ELAST3  consist  of  the  same  components  and  forces,  except 
for  the  rotor  aerodynamic  force  module.  So,  the  difference  among  different 
aerodynamic  algorithms  in  OYSCO  can  be  observed  from  the  results  for  these 
three  models. 

As  can  be  seen  in  Figures  17a  through  17e,  the  results  for  models  ELAST2  and 
RIGIT2  are  in  good  agreement  as  far  as  the  control  positions,  fuselage  pitch 
angle,  and  horsepower  consumption  are  concerned.  Figure  17  also  shows  that 
the  results  for  FRA0,  FRA2,  and  FRA3  are  reasonably  close,  except  that  the 
horsepower  computed  for  FRA3  is  considerably  higher  than  for  FRA2  and  FRAJ0. 
This  may  be  due  to  the  three-dimensional  yaw  flow  correction  and  different 
inflow  algorithm  of  FRA3. 
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FLIGHT  TEST  36A 
AH1G-36A  (S) 


e  16a.  Comparison  Between  Flight  36A  and  OYSCO  Simulation, 
Collective  Control  vs  Advance  Ratio. 
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Figure  16b.  Comparison  Between  Flight  36A  and  DYSCO  Simulation, 
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Figure  16c.  Comparison  Between  Flight  36A  and  OYSCO  Simulation 
Cyclic  Sine  Control  vs  Advance  Ratio. 
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Figure  16d.  Comparison  Between  Flight  36A  and  DYSCO  Simulation, 
Fuselage  Angle  vs  Advance  Ratio. 


Comparison  Between  Flight  36A  and  DYSCO  Simulation, 
Horsepower  vs  Advance  Ratio. 
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Effect  of  Rotor  -  Aerodynamics  Combinations  on  Trim, 
Collective  Control  vs  Advance  Ratio. 
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Figure  17b.  Effect  of  Rotor  -  Aerodynamics  Combinations  on  Trim, 
Cyclic  Cosine  Control  vs  Advance  Ratio. 
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17c.  Effect  of  Rotor  -  Aerodynamics  Combinations  on  Trim. 
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Figure  I7d.  Effect  of  Rotor  -  Aerodynamics  Combinations  on  Trim, 
Fuselage  Pitch  Angle  vs  Advance  Ratio. 
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Horsepower  vs  Advance  Ratio. 
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6.0  LIST  OF  SYMBOLS 


Distance  from  midchord  to  the  elastic  axis 
divided  by  semi  chord;  also  lift  curve  slope 

Blade  cross-sectional  area 

Semi chord 

Rotating  coordinate  system;  also  tip  loss 
coefficient 

Rotating  coordinate  system 

Blade  cross-sectional  Integrals 
Blade  chord 

Modified  thrust  coefficient 

Blade  aerodynamic  lift,  drag,  and  pitching  moment 
coefficients 

Blade  drag  coefficient  In  the  normal  and  radial 
directions,  respectively 

Thrust  coefficient 

Effective  hub  damping  coefficients 


Mass  centroid  offset  from  elastic  axis  (positive 

when  centroid  Is  forward) 

Hinge  offset 

Area  centroid  offset  from  elastic  axis  (positive 

when  centroid  Is  forward) 

flap  and  lag  hinge  offsets 
Young's  modulus 

2  CD* 


E 


v 


E 


T.  J,  S 


Kp>  K^, 

Km’  W  *W 


Effective  In-plane  stiffness 

-  -  EA  ej  cos  $  +  EIZ,  -  (EIZ,  -  EIy,)  sin2  6 

Coupling  stiffness  between  In-plane  and  out-of- 
plane  bending 

»  (EIZ/  -  EIy/)  sin  $  cos  $  -  EA  ej[  sin  0 

Effective  out-of-plane  stiffness 

-  -  EA  e2  sin  9  +  EIy,  +  (EIZ,  -  EIy/)  sin2  $ 

Effective  torsional  stiffness 

-  GJ  +  EBJ  9'Z  +  k\  1*  -  EA  K*  $'Z 

Induced  velocity  distribution  function 
Shear  modulus 

Unit  vectors  associated  with  B  coordinate  system 

Blade  cross  section  moment  of  inertia  from  y'  and 
z'  axes 

Effective  moments  of  inertia  of  hub 

Torsional  rigidity  constant 

Area  radius  of  gyration  of  blade  cross  section 

Blade  flap,  lag,  and  pitch  spring  rate 

Mass  radius  of  gyration  of  blade  cross  section, 
polar,  from  chord,  from  axis  through  c.g.  perpen¬ 
dicular  to  chord 

Effective  stiffness  of  hub 


Generalized  forces  applied  to  blade  and  hub 
translational  degrees  of  freedom 
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m 


N*  ,  My  ,  M, 

xh  \  hi 


u,  v,  w 


U 


UT,  Up,  Up 


Blade  mass  per  unit  length 
Mach  number 
Effective  Mach  number 

Effective  hub  masses 

Bending  moment  about  y'  and  z'  axes 

Generalized  applied  moments  In  blade  rotational 
degrees  of  freedom 

-  m  e  sin  $ 

Warp  term  In  strain  energy  expression 

-  m  e  cos  B 

Generalized  degree  of  freedom 

Distance  along  the  deformed  elastic  axis  « 

Blade  radius;  also  inertia  frame 

Twisting  moment  from  shear  stress 

Time;  also  blade  thickness 

Tension;  also  kinetic  energy 

Transformation  matrix 

-aV 

Twisting  moment  from  longitudinal  stress 

Elastic  displacements  In  x,  y,  z  directions 

Strain  energy;  also  velocity  of  flow  perpendic¬ 
ular  to  airfoil  leading  edge 

Tangential,  vertical,  radial  velocity  components 
of  the  air  relative  to  the  blade 

Total  velocity  of  a  blade  element 
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u 

sound 

V  v 


X,  y,  2 
x',  y\  z' 
xr  yp  2j 
X,  Y,  Z 


xh 


Velocity  of  the  origin  of  frame  B 
Sound  velocity 

Components  of  total  velocity  of  a  blade  element 
in  B  frame 

Axial  force  in  the  x'  direction 

Rotating,  undeformed  blade  coordinate  system 

Blade  fixed  coordinate  system 

Coordinate  of  a  point  in  the  deformed  blade 

Inertia  coordinate  system 

Hub  extent  divided  by  rotor  radius 


n 


Hub  translational  degrees  of  freedom 


Xu»  Xy,  Xv,,  X^,  X^, 

V 


V  v-  V  V 


Concise  forms  for  writing  the  kinetic  energy 
expressions 


z1*  Generalized  coordinates,  amplitudes  of  ith 

J  *  in-plane,  jth  out-of-plane,  and  kth  torsion  modes 

in  Raleigh-Rite/Galerkin  method  (functions  of 
time  only) 


zj-*k 


Modal  functions  used  in  Raleigh-Rite/Galerkin 
method  (functions  of  x  only) 


Y 


u 


/  » 


Concise  forms  for  writing  the  potential  energy 
expressions 
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a 

Angle  of  attack 

“mod 

Modified  angle  of  attack 

aRD 

Reference  angle  of  attack  for  drag 

aRL 

Reference  angle  of  attack  for  lift 

V  v  az 

Hub  rotational  degrees  of  freedom 

^Pc 

Precone  angle 

«(  ) 

Variation  of  (  ) 

sw 

Virtual  work  of  the  external  forces 

AC^ ,  AC^ 

Increment  in  lift  and  moment  coefficient  due  to 
unsteady  aerodynamic  effect 

AE 

•  EIZ,  -  Ely,  -  ej;  EA 

AK 

’  ^ 

€ 

Small  parameter 

f,  V 

Principal  axes  of  local  cross  section  of  blade 

f 

* 

T) 

■  *  ‘  ir 

$ 

Built-in  twist 

\ 

Harp  function;  also  Inflow  ratio 

A 

Yaw  flow  angle 

M 

Advance  ratio 

V 

Induced  velocity  across  the  rotor  disk 

V 

Averaged  Induced  velocity 
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V* 

-* 

p 

°xx'  axf’  °XTI 
T 


*A 

* 

*s 

3B 

Q 

Qx »  *  ^2 

(  )' 
r> 


Averaged  Induced  velocity  with  ground  effect 

Position  vector  of  a  deformed  blade  element  in 
the  B  frame 

Tensor  stress  components 

Centrifugal  tension  integral  [R  mx  dx 

Jx 

Elastic  twist  about  elastic  axis 
Inflow  angle 

Blade  azimuth 

Shaft  speed  perturbation  coordinate 

Rotational  velocity  vector  of  B  frame 
Blade  rotational  speed 

Components  of  hub  angular  velocity  in  B  system 

■i j <  > 
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C81  INPUT  PARAMETERS  FOR  FLIGHT  35A  (CLEAN  WING,  8320-LB  GROSS  WEIGHT,  AFT  C.6 
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